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m.coxjim j $iznm-t&^Mmx&& a^m-^**^ 

m^mtz t <r>£MzmT>mmm&>w.Mm<7)&mT' 
lorn. 

^mn&ftZb<7)imx'M& z t osmm; ttzx^ 

'J y KKIBW*WSil^Wre'ft4i t Srlfti&k-^ 

imbx^m^nm?-? t coft-^^rt <d££#a 

*«#J«^rk<0»CW*i b fcj: OSHttfcLfcA* 

-7 U >y K!E«lcO*IH«^at , &fc*»tft#A:fie 
^<nWr—?\zn LXJ-W LfctiWFMttfc tK^fl^ 

*mw.<rmm>xi}m^<?mT-9b. ztiiznm-z 

?t&?) A^m^six- * *» fesawi * jh v>t as 
Di-mm b ztuznm? a ^oa^nx- * t 
«»r b <vm-izm^wm<mt5 XTm.ftmm&»mF 
<rm*it»0>imvn*-fm Ltzwn 2 mnvm^z* 

<fcfffl*£T-ifc£-r & £ k *»«fc**fefii*lStt?lttfr 

a. 

m^bAjim^nmT-fbnm^mftZbnmz&A 
jim^&ftz b ojaarcw* z t izx 0 aasft a* 
^isgtfj&fr*^ « mmtxijm^^nztui- & jl 

t xxmnnm?- 9 b <rm^mtzb <7>mzm 

mLtzmz&m^f&ftz'b <r>%&xm& £ k tc ± y m& 



itLx^hm^\a,-hm^^zn»h^-9 y f 
cw*3i6 3 mmww<mmi. m-mm^x* 

ii^- b X-h^mnWr—? b (om^fcttZb <7>H£#A 

ximwtftzb<r>'m.xmhzb iz* *)mmitLt:xtt 
gssmtfcwb&to&mtmxjjm^mmz&v&z- 

-9 y -y F8ESt*>gS 1 ^JTiSi&^W&k . *MJffl<7>A 

x^mn<rm=f-9 1 *>m^j*#r t onsris 

■fbLT^ss^baj^fi^stcfcft&j--^ y -y h 
S§Stc7)m2^#illIiS^M|gtk. *t*»t*r>*ASje*o 

mm;. 

^^So«^OA7J«^*f-*-^ b » 

s<oai^«#^«-ri>^a!HiT'fti»A^«-^. 

fs^Wftt^M^itJ^T, A^ff#tm^«^Mf^ 

&zb z&mb^z&teimmi'mm. 
[»*«8 3 mmn-ttwM&u. *mm<?)XJim 
^bxjjm^<7)m7-?b<?)m^f&ftz:b<Dmt:&xj) 
m^&ftzbnjmx®\zz b tci <omm:Ltzxum 
fmfri&w^&to&mfeikxjjm^mstztsvzj-'- 
9 y -y KKgtomi^^wsrr"*) s i t zmm k-rsw 

^k A7jft^o^- ^ t <oe#«»-r k omz&x-ft 
m^fcftzb<vmxmz> z b x oa«fl;LfcA*s 

^ y y K8ESIi0#P^Mitk . t»#A*fi* 

cout 1 -^ izmlx^ Ltzta^mmbiHiim^com 
?-?b<vmKmm*t>9m-t&7-'-?mmmm>t>m 
f&zn&zb zftwib -tm 7 tctesco-fee^ftt 

im&m 1 0 ] mies^wtt^^ia«. ^shi^a* 
m^bx^m^my 1 -? bcvm^fcttzb nmt&x 
■hm^&.ttzbn'm.xmhzb t*«t y«fs^L7tA^ 
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m&m 1 1 1 lTE**ftW>fl«tt. TOW*>A* 

-?vv vvmnm i commm-mmt . ^mmcoxji 
wco&tttm^xm&z-hm Ltz&iim&^tfLMzym 
Lx&fr&MWtiHJjmn&mztm&zL-? u hue 
T-fizttL-zi-MLtiiiijjj'Wm.ttiijjm^nmT- 
ti&zk tmik-r&m&m i \zwm®E8ms?m 

[00 0 1] 

jr^t. wmm*ML, wMm&zmmmwi.t.zn. 

BW=ffli.v&. /yy^-, f^7M* 

[0002] 

htiwmmmz tt:&. ~r v yf-nxjjxh&c m 
YK&mscDm^iz&m-t&mmwmz'ff&K rvy 
?-x'7<) > httiji-t at j*£mi£t& . 

<0*&£\ **^ih-o;iJ;ftT&&3SSK^Lrt:RGB 

^JiOfe fc X* -v OifcfccORG Bfe^EI^fei: <?)W 
flk A#«0fi{2^#tt£IBiKfBfflL 



[0003] #ll^ffi«3^c:fcv^Tffl^& WlflifSc 

gmcoitncoAjjm^ztt-t&iiiJim^zTMLt: o , 

ffi^aj^ft^*ht&A*«-t$:^S!lU^ 
k&S. 

[0004] ^-AattiSa^feeg^^Krfei 
LTJi. Aa^ScO^aW^ttSriESUT^^/Wk 
U^JI*T/P£J^T^83-$-Sk<Ok. ^g^Affi^ 

[0005] #^F5 - 1830 5-§-&?8T 

l±. A*SST'2bS*7-Xdr^i— t;iJV^T{±3X3 
Vb'J iiijjmWX'h&TV yf-lZ&^TU 

mmmmx$>^> l. ^/u^m^xttzk vxtm 
[0006] zttizMLx. 3yei-W5-77f 

St^^L, 1»SOA^«^SS7-'-^t-e^»i5-r 
h tti^fi^O^r- ^ cox- *b WROflHRfcSBe 

#{c^tf-sA^m^&^}»j-rsi.^T"$)&„ u>>u % 
?—9<rm ^k^mzmt^izmumwrnm^ 

[0007] 4fc, #^7-8 7 34 7#^fS. t>J: 

, m?mz&£:m. vol. 35, no. 2, 1 

996*. P. 125-129 izii. «|gtc0A^fi-^O 

fi^i=str z &*m^ z^-w Ltz <o . imoixjim^iz 
tt^&Atim^-zfM-r&iTmtfmmztix^h. z<r> 
■mx-a. tmizemmwi-z'&kfr, zmx^m 
mkmi<. m?-?v>&£.isi<7>j>\-z9mLx^m~?z> 

[0008] ist. »if s P2-2 2 6 87 0-5§-£r?g-C 

mwAjim^mT-'-fkztiizium't&mjjm 
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«IJ»^3««3EBI*l8»tf5!rctSVfttt, fiES^aj^J 

£tz. ZttmiiAJjklhJwttTti)* 
niKTtn^lZ&th^mX'b^X . AyWYMCK 
d^4}X7C-Cffi^A { feSriS^-r3<J:7CC0^^Ji^ 

[ 0 0 0 9 ] a , XSn$Rff5 4 7 1 3 2 4^HJ§«« 

k^xta^m^t^xijmm^mt&^tizi*) , r 

-fflco^/t/ y ? T -y r^-^c^mSr^S 

2-22687 O^^gfctJfflSia^S* 1 *^^. 
[0010] £tz. ftffl¥2 - 2 8 9 3 6 7^&fg£ 

\±. m&nx^m^nmT-? tzmznm-h&Jim 
wmr-fw-fm**, mtwkmmztihtm 

Ic^ttTfc 0 . ±3£<9!|f H¥ 2-22687 0^$Ri: 
R**>IBHKf»*4j&\ S>&^ ii. &tt£.mizftW\-?hk 

[00 11] 

zm-mzmw. i xtemmm^mt h-n 
-t%*>*>. &n<ri\tim^zn?h&?imm:*m 
-th^m. tixif. &n<vauim^izMi't>A?>m^£ 

m^bAjjm^nnvn-gtizTmt&iTmkLx . w 

k£BMki-&i><7)X'$>&. 

i . mm^mm ummm) fctt»L*v^-r*4 
i: lt , xj*-is>7mm&ft^>ttX'bz> z t . 

3. A^ft^-^TC^MiJGccO^Tt. 3mt<r)t%r& 

tmmzwiz.$,zb. 

4. : FHtftKjW«V*£4:. 

5. ^}RJ<ic0iiM14*«1«Cil$^-5^t. 



6. 

7. faw-a-?**7 bV-^fcfflv^^lcibK 
[00 12] 

oft, j m.^x-hWA<n i m—9h. zfiizMmthft 

fcr-fixjf. '<TW0f&ftZmfet&Zk1) { X'Z&. &tz. 

n*ttwimk Lxwrnammpt thizmicom 

* i: com^&frzk <nmHzm*-imcr>mnTm*ttv 
Lfzmto 2 mfatfm/nz% h <t a tc . fr*i*>j£# t 

[0013] *4 Wi. ifiMf^ffll. A^3«^ 
Sr^SMi: *5t(±ai^ff#i:-a5^A^ft-f-Sr«^ 

mm. wwAum^z^mmtL. xnm^wmrr- 
^^^mmm^x^tzta^mmk wznm-t 

iz%& xoiz, 'nmn&ftk m*.ttw>im.k ^mwz 
m<m{mxikiz-r&zktfx$2,. zcom&. mx-n 
ttcomnk txwf-mmmfk t ^zm\<r>mtm^> 
imzm^xtx^. 

[0014] ttewrtwHfei* . mwHMN&b t 

X. h^frtfr&AXm^mT-flzKLX^WLt: 

ta^mmkta^im^m^-^ t <vmmm>t>m(- 

[00 15] 

d ^ x-it, Aiim^tf 3 fc«4 jjscrajsifi-w 
3dttv®&zmzjR*), OTt^-t ( 1 ) 5W>>£ 
( 7 ) ic*jav^»w-*-&*», Atitii&ttdsatttz 
ftt>izmfe2n&*>vx'%< . fa^7ci&Ttiii«-cft 
■s.. < 1 ) *k ( 6 > A^fi^3<):7uT-ai*fi 

^3^7C<7)^-{C*fJiELTt5 0. (2) 35k (7)* 
fi. A^ft-^*^ iJ^ccO^C^JtE LT±3 0 , ( 3 ) ^ 

[00 16] 
[^Cl] 
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{Bbs i = 1 ~ n 

Ei = E(^((yi,-yH) 2 + (^,-!fe0 2 + (j/3 i -y3 i ) 2 )) 

+ E(^(yL-2te.) 2 ) 



(3) 
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[0017] *r . ( 1 ) *t ( 2 > *l±. nifi<?)A^ 
ft-^-iO^-f— ^ { x u . x 2i . x 3i ) (xji, x 

2i . x 3i , x 4i ) , i = l~ni;. -fO^iM 

<yii' . y2i' » y3i' ) . i = i~ni:£, SeSPS 



[ 0 0 1 8 ] #0 ( 3 ) sStf. IS^A^ft-^^Hx- 
( 1 ) sSifcti: ( 2 ) sM^^MSrfflVvC;^ 

y«i' . y 3i ' ) . i = i~n«i. 

(2) 5er^t^nffl^A*fi#«0Slr-^ (x, it x 
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2i p x 3i ) ifckt ( x lii x 2i > x 3i > x 4i) > 1—1 

(yii, y 2 i. y 3 i) . i = l-~nl±. (x li( x 2i , 
x 3i ) ifeti (xu, x 2i , x 3i , x 4i ) , i = l~n 

[0019] ( 6 ) 5fti. AJ)flW3<X7c:T-ai*fi^ 
Jtrntr^LXH*). mu, mi,. • • • (2 ( 1 ) sSkPI 
23> x 3j ) (6) ^(^A-r&^t^«fc 

"9, mmz^mmx-to&itijim^ (y^, y 2J . y 3J ) 

4. *=HWW*i«jf&f - (y u . y 
20. ysj) T"S)?»^<2, BMRC (6) ^S:^»<i 
ktciO, TOWTObtAAfi* (x u . x 2d , x 3o ) 

[0020] ( 7 ) 5N±. A^fiW4<X7C-C-ai*«^ 

„. m a2 , • • •\i(2)3tbntm\<r>&®.ftX*fo 
6. WL z fWmt i X-h<& 1 % (x u . x 2j , x 3j> x 4j )T 
£>6*ll=r<;L (7) ^fcftA-tSitfcJ: 1 ?, SWCT 
«trc**ffi*fS# (yu. yjj. y 3j ) SrsftftSifc 
Wf-MmsWiJ*©^ ( yi d , y 2d , y 3 j ) "C 
&&*§-£<i, (7)^^i¥(c»<ifc{±-C-#^v^coT\ 
A^«#0-S5, WxJf (x 4j ) WP3SHik LTt!t5£ 
L, ( 7 ) 5S£i£Kfi?lvC, ?SO^afflT$)i»A^fi 
#f»*A*«-»f, flRlf (x u , x 2j , x 3;j ) 

[ 0 0 2 1 ] CicT) J: d tc, A^m#0&?T#4&7E?)*§ 
3^<o^kla|^fc:8id.rka<T*#S„ 

mmwmi/^»x^htzih^mmtiim< , as 

^JcOjlx-^^Aaj^feS^O^Sr^^'-LT fcvfr 
< T t> A*^feS^±^te^T«ffl«rSfC'35 0 . A 

a. 

[0022] *W2*£kOWilS^O:t^£ffl^S 
£ k lz£ ~>X . «*ftft<0«MRWi j . i = 1 ~n swea 
-jTV^^tcEj £g/jNk^£«fcd&^-CfiW# 
em u , mi 2 , • • ■ fcfc' 



fcii^mrc* & A^ffltcte#-r^^<i-««tc 
i&b&fts fctfreupfirvwe. (3)^Ej s-s/hk 

-f h fc jftftoTT-. fi^+ttO^fcWij , i = 1 ~ 
n fcfi^JCOJ&ttmn , m, 2 , • • • k ^SlHfiOSj&ffiS: 

[0023]^, Ej ^ft/hk-rakv^^kii. 

( 3 ) XZgBLti ( 4 ) 3K.frt>bfrZ J: ? fc. HJ^ft 
<nx\ '&«<nmnmftlzftMLtzm.*W3LX'Wt>hb'?& 

t^ozkkmmx'hz. &o(5)5W±. mmx-n 
m^nm^-ffr^nmzm^xigtbf.ziiifs^'mmk * 
iiizttm-tz>wm<oiiijimT : -?<miim-tfittz'k fc 
%.$tmLfcmizttLxnA-ttif2tiK$m<?) 2m\ 

Ej ' ZijkLX^h. *m\coljmii. (5)3tlZmL 
tz£o%Ej ' *&Hzi-Z%&lzi>^<m&tZ7m 

wirmmnxjim^ t a*<i^-<oht-? k <7>m^& 

mrm?- ? k emmmzk <n^mfn^^m 
i ^xm&ztofctti&xmRcv&MzimiLti&tzjm 
itLxxfr&Mteit&jjimgmizinf&jL-? u v f 

-ffl«AtBH^aco-feMa<?)^ISSr^^ i k 

[00 2 5] A^fi-§-A J 3»5c47t{24^7C-t' 
ai^^S^TCCO^SrPl^i: 0, (8) 

(12) ( 8 ) 5tk (1 

0) ^{i, A*«^* i 3^7CCO«^»J5eLTJ3 i y. 
( 9 ) (11) 5«li. A^mW4»5cT-ai^lfi^ 
3^7n^^-^jS LTfe 0 . ( 1 2 ) i^ti. 

^{znmLx^z* 

[0026] 
[*4] 



4 
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W U j = Fij(((*u ~ *i,)/*io) 2 + ((«2i - x 2 j)/x i0 ) 2 

+((* 3i - x 3j )/x 30 ) 2 ) (8) 

Witt = ^u(((^l.-^li)Mo) 2 + (( a:: 2,-^2i)/x 2 o) 2 

+((a?3,- - x 3j )/x 30 ) 2 + ((x 4i - x 4j )/x M ) 2 ) (9) 

W 2ii = G 0 (((m n (*« - x^)) 2 + (mi 2 (x 2 i - s 2 ,)) 2 
+(m 13 (x 3i - x 3i )) 2 )/(y 10 ) 2 
+((m 21 (x u - xiy)) 2 + (m 22 (x 2 < - x 2j )) 2 
+(m 23 (z 3l - x 3 j)f)f{y 20 f 
+((m S i(x, i - Xij)) 2 + (m 32 (x 2i - x 2 ,)) 2 
+(m 33 (x 3i - x 3i )) 2 )/(j/ 30 ) 2 ) (10) 

W 2 ij = Gij(((Tn n (xu - x xj )) 2 + (m l2 (x 2 i - x 2j )f 

+(mi 3 (a; 3i - x 3j )) 2 + (m 14 (x 4i - x 4j )) 2 )/(y w ) 2 
+((m 2 i(x lt - - xtj)) 2 + (m 22 (x 2 < - x 2 ,)) 2 
+(m 23 (ic 3i - x 3i )) 2 + (m 24 (x 4i - X4 j )) 2 )/{y 2 af 
+((m 3 i [xu - x xj )) 2 + (m 32 (x 2i - x 2 ,)) 2 
+(m 33 (x 3i - x 3i )) 2 + (m 34 (x 4i - x 4i )) 2 )/( y3 o) 2 ) (11) 



W i2ij = H(W lih W 2ij ) 

[0027] Sfc-T. ( 8 ) ( 9 > *I4, S*tttt<0 

(X U , Xtj, X 3 j) ifcJi, (X ld , Xjj, X 3 j, X 

4 j ) * { !S-f-fflI<l^^(±^}fflfficoA*fi#. ( x liP 

X 2 i, X 3 i x 2i » x 3i» x 4i 

ft-f-^-f— ^. (Xio, X 20 . X 30 ) ttzlt ( Xio. 

x 20> x 30 . x 40 ) #A#fl^£ia<0^,M:<05gfa. F 
ijjWWWfcHWft-C**. «itf (8) tfX-ii. M^m 

mtKtej'MmcDAJim^ (x u , x 2J , * 3j ) tx-x 

m^<F>Wr—? (x,i, x 2i , x 3i ) iOff^frri; 
C 7 )^ ( Xu-Xij ) , (x 2i -x 2 j) , (x 3i -x 3 j) 
ZthZMffltcVjm (x„. x 20 ) x 30 ) X'M 
t&itt. (xjj-xjj ) /x 10 . (x 2i -x 2;j )/ 
x 20 . (x 3i -x 3 j ) /x 30 J:^"S. <I*i.£<023fcftl£r 



(12) 

[0028] <jrcc. ( i o ) (in ajti. amni 
fi^fcti^PJffi^A^fi-^fc A^fl^Hx-^ t com 

S-^LTt^. iit, (x M , x 2j , x 3 j)ife(i 

(Xij, X 2 j, X 3 j, X 4 j 

AlliE^r. (Xu, x zi> X 3i ) *^{i(Xii. x 2i , 
x 3ii x 4i ) 1fiXK^<F>Wr—9 . m u , m 12 , • • 
• im^'irn<r>\m . ( y i o . y 2 o . y » 0 > tfaWrtl 

*mw* Mttmmcoxjjm^ t x^m^^T-^ 
t vm^f&ftzt con t w<i<v&ftz®VT fi#f£#r 

fc^23HU^c»S:fcS^i;^«i*UTfe , 5. 23H2: 



(8) 
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ft* k &®&tziz. &m<r>&mz X otlia>«t«l 
<X^\±A^<^-yX L4 5 (10) 

&j:xs (id 5H<v£oiz2m-i>z.nz£'>x&m<D 
tmmzm&k txtm-tz z t *>t# s . 
[00291 mti* do) 5^2. m?$m&M±^m 
m^xtim^ ( x ! j , x 2j , x 3;j ) tAMf<o^f- 

? ( x u , x 2i , x 3i ) tOfi^^i'tcoH ( x u - 
X]j) , (x 2i -x 2;j ) , (x 3i -x 3 j) £5£#>6„ ^ 

fucfi^gNts z i: t j: o m^aga^js^feisaiiT'* 

fcftfcWfcfc*. (iiiiilxu-xu))' 
+■ (m 12 (x 2i -x 2 j) ) 2 + (m 13 (x 3i -x 3;j ) ) 
2 . (m 21 (x,i-X]j) ) 2 + (m 22 (x 2i — 
x 2j ) ) 2 + (m 23 ( x 3i -x 3j ) ) 2 , (m 31 ( x,i 
-xjj) ) 2 + (m 32 (x 2i -x 2j ) ) 2 + (m 33 (x 
3i-x 3j ) ) 2 £*l^t«fftfliOj@R 

(y I0 . y 20 . y 3 <>) ZRI^XZtiZ'ix (y 10 ) 2 . 
(x 20 ) 2 , (x 30 ) txm%LX%mtL. MM.LX 

mmm'jm&GijW^t-fkLx^z. 
[0 0 30] mz. ( 1 2) awi. ±.&c?>2mmmts, 
'j?mL*iii8.Lxm*wi<?>{m&fti&'ti> z t 

mk&m(7)£dzmmzm^zzktfx'Z&. 

[003 1] Z<7)£otzLX (12) ^X-mtb^tifzW 
itutf (3) 3«<0«*ft»t<0flaKWi j T'*i. -1*0*^ 
ttW>ffiKWi j £ffllvC&/h2ft£(::J: ( 3 ) *C 
ij^Ej ^S/ht^fif^JcO^mn, m 12 , • • • £ 
3<tf>&. -etT. ^fcffWOjR^m! L , m 12 , • • • 

frt>%h'fim*m^x "f-mm&smi-z ♦ $ 

jru mm<o^mm*nz> z t *>t# § . 
[ o o 3 2 ] si i « s *min^^^fwm<m 

S 1 1 cOAft^Hx-^ttWXS 

o^O. BuifrO ( 1 ) - ( 3 ) *T»<fcn« 
OA^ft^afr— * (x,i, x 2i , x 3i ) , i = l~ 

n k s *iuz*tm-r& nm^ta^m^mr-^ 

(yii. ya- y3i) . i = l~nS:J8ft^"&. ^eofl 



£ 805t t T *> h m * CO * 5 --"C -y f- £ S^ii * ■* T tti fi 

* b* k, mu<r>RGBmk<Dm7-**f*mm-hz 

ZfflfethZklZX*). A^ORGBffii:. ftJ}<V&m 
Ss L* a' b» i:c0^f-^^£fflSt-S£ 

^ Y M C K %<7M.* COm&tX'* 7 — /t >y f-£ tfj;fj 
L s ^cOfe£^-r&£fctCJ:0. A*)OYMCK% 

k. &jjcr>&&mm. mm. l* a* t>« kcom?- 
fnzmm.'th z k izttm-rz . 
[0033] s 1 2<Dm^mmmmrmx. mm 
<n ( 6 ) sJtfj ( 7 ) SCT^/ilR^jBMI, o* d . 
S^Mffl^AAfi^co^ii (x 1Jt x 2j , x 8j )*3t 
fi s (x,j, x 2d , x 3j , x 4j ) £ s IS^lSUfi** 
as*S#^«^i. (y u , y ti . y 3i ) . 4fc«. 
(yu. y 2 j. y 3 j) k <x 4J ) £. ^SSMrSWfftffl* 

[0034]»:. si 3»a^#(t«of^co*JWffiS: 

mfct zmizx* h mr&miz-t h z t 
#>htz#>[zmmx'h& . «i«r, ^suwA^ft^-o 

J^-^li. ( 8 ) aSfc ( 9 ) ^tc^L^life^ mmXhIMft 

k A^m^^Hf-'- ^ coft^js^r k co^mmt ttz 

fflv^s. z<r)fmxiiffm<7)jsiftwm5£ztLx\^\.v?> 
x\ do)*, (in *tc^ Ltzmim&m®G i j 
i±m^-r. «*w lid co^T«<at-ts. 
[0035] m-mm#iiiJim^<&&isi±. tirtgmx- 

^^L^te^-suffiaj^fi^coiiT-^ kcom^f&ftz' 

G i;i c7)a^W 21d '$rfflv^. A^AMiXTLW^- 
{i. (13) st(c* ltzm*w iU 'k. (14) 

u/iA^^ia^ 1 ^<nmtfmt*m®m.vK*3}m 
^QmzMf&jmn&im&m&G^zj; stw 
2i j • k £, (15) ^.x^-rmnHX'^Ltzm^N 

[0036] 
[St5] 



(9) 



10-262157 



= Gij(((yu - j/ii)/2/io) 2 + ((to - 2/2i)/2/2o) 2 

+((?/3i - J/3i)/l/30) 2 ) 



(13) 



(14) 



[ 0 0 3 7 ] mz. SI 4(7)tf?'Jc7)J«^ajXS-C. 

x-^fi^ftftSft-fci-? 'J «y h'3g8i<02 JltWEj 
T, ffM^lft^mn. m 12 , • • - *Vb i-fHW 

[0038] KK, S 1 5cr)^fflffi»atJXSt\ ( 6 ) 
ft&fzte (7) «™ia><b^$H6SrtAfc 

^^jra^r*. *!t^aijffi* { A*m^ (x 1;i , x 2j . 

x 3j ) 3.fzii. (x l3 , x 2j . x 3j> x 4j ) 
(y u , y 2j . ysj) S^affliVaj^fi^ (yij. 
y 2 j. ysg) ttziz (yu. y 2 j. y3o) ^ (* 4j ) ^ 

i^(± (x ld . x 2j , x 3j ) Sr. S14^fr?iJ^^* 
ftX^T*^ff^g-lffcfrm n , m 12 , • • •Sfo 

[ o o 3 9 ] si 6c7)s^#(t<7)f*iSiositaix 

ST, ( 8 ) ^4 ( 9 ) 3^ 8mHffi *fcii^88 

©i*w li3 t v do) it^^ji (in ^ s «ff- 
m^^-r fc o*^tf jflois^s- ffl^rmg^ani* Ltz 

GijOS^ij (12) tt&mmnX'it&L 
L . n ?mzft Lxm^N, 2 i j 5rS«tf»^- 

[ 0 0 4 0 ] SI 7<7)ftP]<D!&.ftniliJMX\ 

( 3 ) &x'fnLtzimcr)AJim^nmT-?frt>?Tm 
m^x^fzm^mmt ztuznm-t&immiijim 
f-rcDmwmztifzj--? y «y pffigico2^WE j 

£«*Cc-f S i: V> d^ftcOt t ®/jN2Si}££ffl^ 
T. *T?<J<7)#J&frm 11 . m 12 , • ■ • $rHJK»:jti-tl>. 

[0041]^:. si 8<r>l>mmi»liiJMX\ 
( 6 ) sSifcli ( 7 ) sSfcfflwt, ttTOM*6TOHI 
tSStttB-ri). JfHMWAAfi^ ( x Id . x 2j . x 
3 j) 3.fzlZ (x^, x 2d , x 3j> x 4j ) <7)i§^4i 



(15) 

(yu. y 2j . y 3j ) mwmtf&Jim^ ( yii . 

y 2 j . ysj ) itzli ( yu , y 2 j , y 3 j ) k ( x 4 d ) o 
(x u , x 2j . x 3j )2:, S 1 7 cofimcOf&ttW 

ntaiMT^tiffm^^m^, m 12 . • ■ • £® 

[0042] mz, S 1 9 ^i0fflcOlR^I^XST\ 

s i a^mmMMtaxnx'^ibtz^mmt < si5» 

T-fflmHtiJM* fctt 1 Hlmr^S 1 8<tfFW«B*fflI 

vmwxmt lxh. mm. 
mmtft^fflmtcomcom-z&L-ftztizxiib. 
mtmtffot>frttbm%.Lx*}^fzf&ftz-:bcDmf&£ o 

/h $ V S »iiR« L T V ^ fc fij^-r &z\b ifiX'% S . S> 

ZlitimLX^&b'Pm-t&Zb >1-T\ JR 

^LTV^SkfJ^L^^^if. S 1 8cO^MffiS*ajX 
gT^0**fe^iMi^S*W^ i Sr^fflffifc LT , S 2 2 
co^i!l^TMSXg(Cjt^, iR^LTV^^tfl^tJt 

isr^if. s2o<r>?mm<ommm%jMizmtt. 
[0043] s 2 o co^mm^mmm^iMii. 
fcvyuwift b wmnnmiLzm*.. § h izwmm 
z'n%k>&hi)^frzmm^&iLUx~bz>« cwxmz 

L^m^iz. im^<m%>v-y{ l zm^£^xoi l zL 
x^&. mm, i coxs^i Ha ^frt- 

ffii: 1 Wto*tfmb<mH&XVimtteyfBMb 2 
m%&1jftlzis7 hLX^h b # (i-ecoiu^^^^y h 

*l^««-«UfJ:^. dd-C. ^ail«t<7)Sft«SrfT^ 
3 t*BELfc«rfetf , S 1 6c0*^#(tO^ifc^W»tii 

jim^m^x-wmitzmrfML. mim*ft%&%^ 

k WeLfc=5r^«, S 2 1 «>TOfiglU*&£X&&S 

[0044] s 2 1 <nj>wmmmmzxmn . ^fflffi** 



(10) mm*? 10-262157 



h &mm<vmm.x~mti ix i & -> x mm l&v 
x't>h. nft&z'ffK-itz^mmzmteLx&z . ma 
sua. j-mmt lEB^snffltioii^fit/h^v 1 . 

[ 0 0 4 5 ] S 2 2«^»ftTP J£lS«, itO 

fi^f&^mwmi&xm-T-rz . 
I o o 4 6 5 act:, #^HJ!«fefE*^tt^M:fr&^2 
cOWtOJBBfcoVYtBWH-*. d<73Sf£2<7)S!Sfc7)762g 
T'(i, ( 3 ) jZconX-ttWf&mVn t LX . ±^<om 

mm* wi^mmnKhm^ t xjim^n?—? t <o 



»Gij(0l5a»t=. ^-^ffiS^SJi Srmv>l>0lJ5:^ 
itLX±Mi2W&)**W 1 is kW 2ij frflgfflLT^Hfl 

m*hLZzmi>-t%><r>x'foh. 

1 0 0 4 7 ] £ £ X\ ttfiWti fcttTOWIOAAfi^ 
coHir-? k co«^fi£#.r'k Sffift L*:A:fi£f8J 

m*mnm<Qm-£m^x®m£tomLK&jjm®cr>i& 

fcfcit&JL-:? >J -y Hffiiffi^mSfljS^BgStGij k«i. ± 

mcom i <^mm<r>mm Ltzi>e>bmtx'$>&<7)x\ 

xi)mnnmT-?izftLxztit>cr ) 2mm<r>m.?>-$:®. 

mLx=pmzti*:}iii}^mmb\iij)m^<7>m¥-?b<7i 

6 ) 5fo&'4> (23) 3CS:fflV^-C»i»r4. 
[0048] 
[»6] • 



(11) Stmn 0-2 6 2 1 57 





= y'n - yn 


(16) 




= y\i-yii 


(17) 




= y'xz-yis, 


(18) 



±Vi« = (EWi'* - *»*))) / (E( v *>) 



fiu mam 



(19) 



(20) 



(21) 



W 3i = Ji(/ii (Ami, Ay iB i), /.-2(Ay i2 , Ay, s2 ), ^(Ay^, Ay iaA )) 



(22) 



[0 04 9] (16)^(18)*li, imBcvm 

y i3 " > tli&imT-? (yn, y i2 . y i3 ) 
SSr^LTfc 1 ^ (19) (2 1 ) stti. i#B 

[ 0 0 5 0 ] ( 2 2 ) (16)K (2 1) 

y i2 . Ay i3 , i#@<O^T-^^S:^<ll'r-^MO 
tt(t¥%^l^#fik*c9ilO*#£Ay ial , Ay 



(23) 

[ 0 0 5 1 ) ( 2 3 ) iSfi. 2o<7)*li^KiatF 
[ 0 0 5 2 3 H2li, ^^WfeffiitWttWWrttOflt 

m<mmtmix'foz> . s 1 i<7)Aai^iir-^*fii(i 



(12) 



0-2 62 1 57 



4 ofc < R***>?, SMHtt**** . S 2 3 Of 
\, vta*ftftOfll» £f WTT & . flffl 

L$rv**&fcU:s i 2N&&, ±^(7)^ i <r>mtcr>mm 
RiwwiiiNtatf. 

[ o o 5 3 1 s 2 4 (VT-*%im&M*m$i.miim 

3fc*3*Vtvv&. J2TF, a3S-ffl^T^-TS. S3 

i coa^«ffl^si^-c\ H2<osi icr)Aai*^x- 
**ras«xe?*» Ltz^x nx-tiWf-* zm-wm 
tvxm%.L. S32 ot-ca^Ht-'- * £*ti-& aj^i 

[ 0 0 5 4 ] mz, S3 2coa*tttt<9^cO*#JJfetfc 
?£XgT\ ( 3 ) nfflcOfi^^ttiO^Wij 

(Tmrnm^ta-t a . «miafcfc. ( 8 ) L^ts^ 

«-9«IBIKiJ»t* *.-9 'J >v Hra»0*i«S5WJ»Fij 
[ 0 0 5 5 ] ifclc, S3 3c?5tT?iJ«7)^^:ajX^T. 

T. ffiaj<7>£-J&»m u , m 12 , • • • SrlAi: 4-f»dJ-r 
ZZX'<?>mfrWJ<7)fc$StW u li. S3 2coS^#ft 

10 0 5 61 WZ. S3 ACD^mSM&lMX'. (6) 

js-^asffl^A^ft^- (x u , x 2j , x 3j )t*sOT, 

(yu, y 2j . y 3j ) S3 3^ffW(7)jft4Mjtajxe 
T-*tf>5t^iJO#R£^m, , , m 12 , • • ■ m^XHtib 

[ o o 5 7 ] mz. S3 5cos*tt(t?)^(7)B*an: 
st, (8) ^iZTp-m^mmttc^mmnoxijm^ 

fcA*fflia^«4h&»fe*ft*ai«bA*«^aaifc:*j 

VhP--7 >J -v h'gEfitO^^^F i;i c7)fi^W Uj 
h , (10) ^fcS^H^ttifcll^HWt^A*^ 

mfrzm^xsm&tmLfcftxmsnmMz&m 

(12) ^^■T^raifcH'r'-^fiKL? ,k v:*^w 12ij $-. 

n«««f-^»K« tx nmMmt h . 

C 0 0 5 8 3 mz , S 3 6 (^toJfi#?|*ffiI« , C. 
(3) 5CT*L^1ilS<7)A*)S-f-^7 ? -^*^^iJi& 



T. fiMOttKhn,,. m 12 . • • • ^WSJUft-fS. 
[ 0 0 5 9 3 <XtC, S3 7<7) : fMtmniiiXWX\ 

(yu. y 2j . y 3 o) s3 6<o?T*ijojfcfrgtisxg 
[ o o 6 o ] mz. S3 8co^aa<a<?)jR^*ij^xgT. 

S 3 7 <^«9^J^ffiX|it^fc^Wffi b s S 3 4 
c0^11ffl*JiiXS4?t(i 1 0HffOS 3 7^PJfflWEHi 

xmx^^^mmtumLx. ^smwoGRui:** 

WSfcStttffci: LTJi. 0U<£. 

HSajXgT^^wt^SlffiSraiKW^aMtiiLT, S 
4 i^»T«^xS^jl^. «JKt*CV»«rv»kT!je 

LJt^r fetr. s3 9<«««oiffwweijefc«tr. 
[00613 S3 9<r> : fmmconnmmmxmi. f-m 

'tth^&t^t^mfet&JMX'foh. z\<mmizx 

^u^-ir. mmcotfs.f^—Titzm^Kc^xdizLx 
v^-s. £cr)is*<iinig<7)^i. i^fr-c-^ 

*w<7y?Mm<nmmimzb&m*, mt^mz^y 

WjhUfalzis? h LX^Z, t^li^St^^^yh 

vxv<r>^\&m\nzintc7 bmiZL. wmzmatzb 

fc^fcflJ&rMHr.tv*. ££T\ ^aiffliOHti-SiSrtT^ 

3 fcwsLfcssr^tr. s 3 5 <r>m^mL<mn^,xn^ 
mA,xmm**mmL* mt&&ft%to%^bm. 
Ltz% s 4 o <m-mmmim^i.mzmts . 

[00623 S 4 0 tfO^aiMiSi^^XgJi; . "HHflb^ 
«> &fiHlWJ51fflT« LT L 4 X JRsg L=5: ^ ^-^tc , 

mh*nt> L\,m**mmrfmmt Lxm%.-t&xm 
mmb^ftbmM^^mmmi-thKbizx 



(13) 



0-2 6 2 1 5 7 



[0063]^S41 ^MRTHSISWU iT» 

[ 0 0 6 4 ] JJW) S 4 2<0^-*»JfflH»«**aSie 
(16)^(22) 55fcyj?Lfcf-^lil*BWt 
Ji CLt*»ot, T-^WKWfcfi^WsitJMJ-r 

*<*>m0rjz § 3 *» 6 . WfrtfH fc*KWf £>H«9A * $ 36* 
*£<3:£15k'*£&<I£. JWr**i6|«r4>«a<0^:# 

[0065]^S43 ^7*-^»«M««»J|tffl»T 

araatc «fc *K**asa**T LfcjW5a»*ws&*-*ie 

T\ 5fo»T<0*aii. S4 2*)^- JKti 

sp^rka o>g-r. &Tw*^Ji, 02 

cos 1 20*mB»MtDBKatr. 
[0066] E]2^S 1 2*6S2 2£?<0?flDj|l 



[0068] SI 6<^SA#lt<0fldK<0f^mxa-r. 
( 8 ) a$4fctt ( 9 ) *tC* , «ff^4fcl±^ai«0 

ii^Wnj k . (10) (li) iWc^-TS 

^»ffi4fcti^a!ffi*>A*ft#k A*ffi^>#7*-* k 

raiScGijlcJ:Sfi^W 2ij k. (22) Mz^Ltzf- 
^»KW»Ji t=i*ft*w 8i fr. (2 3)£*>*jftN 

[0069] JjUitfDS 1 3 k S 1 6Ja*Hi, T-*»Jg 
0. s 1 2HS2 24-W)fl»a*at:i:l:J: 



at, T-?«M&ffi^£ffio^^ tt*><fn*&«> 
1 3ff>m/f-iHf<omL<owsim»tiixmx\ ( 3 > ss-c 
£*>*j»nifcL witf . tfc*ag»<A*«*«>*£tt, 

< 8 > sfcifcfi ( 9 ) 3Sfc*LfcBTOHtA*f3#i: A 
^fi-^^^T-^ k ottffllftr k co^S-SI&ft L/c A 

u y Hm^nifeHnR Fi d cos^w, i j k . 

^-^HJWWUteWjifcOTfcii?*- ( 24 ) stT^Jfc 
UfcfcWfcJH^*. »^Mfi* s tti;>J<i^-coi§-£{i. A* 

^3<j:7c<7)^-{i, (13) ^ic:* Ltim^wmxjim 

#k (ttflB**)*-?- * k W«#JS*rk <OHS:*t&fli 

k , ?-m$m§m*-w3i t mnt^ki- ( 2 

5) rtTfeJftLfc&ttfcfflVi*. A**M#tcco*§£- 
li. (13) W Ji;j ' k . (14) alSfcwLfcA 

1 o«o^Bg^*tt*a^kUfcA^«#^ia 

(=fe(t4^»*iBftJdWRFid Ki SS^Wu j ' k , 
7-'-*««Wifc£*W 3i £. «Tfc*t (26) sCT^ 

[0067] 
[*7] 

(24) 



(25) 



Wai) (26) 

[ 0 0 7 0 ] mz. *9t&ttmmwtrfmrmPM3 

Tli. ( 3 ) JtOfiRRWii k LT .SB 1 

*<^B«-CJBv^*mMS^)A**-f-k Ajim^co^? 

&Z k tci DIML^»»«ftS» 

[0071] zzr\ m-mmx-^m^k xh<m<r> 
sij^^ k fc a "3a»fl:UfcA*auaia»**3&»«.*»4 

^-MiicF i;j kJi, J^^mi^H^^f^L^ 
(8) (9) *kntT*«. *HMT?*« 

tmv>&+im^iztt&^&^mm'Qhz>Ktim^*i&#> 



(14) 



0-2 62 1 5 7 



im-r 4^rc eLtsmnvti-wirmi. < 3 > ^nm* 

#(tO#StW i;j tU s (12) *W 12iJ £ffM4tt 
(8) *4fcl± (9) ^(C^L^Wi^fflV^S 
J3l*Wi» Blfciir«r8ci:£<H&fca4. 

t 0 0 7 2 ] H4« s *HBBcofe^#tt^iffl*&<OSI 
3cOlfetec^E»fcirc tJfHHrc*4ffiBtt>A#e 
^tt^4^ffiT&4airtfI^£#a64«1^fctt 
S»f^-#J£^-t:7P-?-+-hT-;&!.. etofi, S 

i (c^-r s 1 1 coA&jjmT-fttmffiTmt^K mm 

tew, Pllif«tS. 

[0073]tf;t;, s 5 2c9«^Wiie-m> » 
^ll!lffl* { A^m^^^=5rcOT\ ±MV (6)5$. 

( 7 > XTm^m-mmn^m^ ( Xlj , x 2j , x 

3 j) i/tJiA^fi-^- ( xu. x 2j , x 3j , x 4 j) 
[ o o 7 4 ] s 5 3<oa;ww*ffcD3:an:s 

T\ ( 3 ) A-CSfrtfc nffl<OS*Wt<OflaRWi j £ *iti 
-T4. S^llffl^A^ff^^T'ftO. S^tttfcDtfi 

ist lt«. ( 8 ) 5t*^tt ( 9 ) ^ifttmn 
oA^fi-^t A^ft^o^f'-^ tofi#^rt con 
vrnx-mz z t tct oaigft l*a 

< , — 3BWC*«> 4 £ fc arc s 4 . 
[ 0 0 7 5 ] S 5 4<DftMe>f8.ftMiiiJLmX\ 

( 3 ) 5CT'*t^1SS)tOA*(i^llT-^*^tf^J$r 

[0076] &tC. S 5 5c^$HI»£tJX^T\ ( 6 ) 
l£&tzl± (7 ) rtfcjflwc. S^SM*>^$Jffi£&ai 
-T4. «T«B*«A*«# (x,j. x 2d , x 3j ) tfzii. 
A^fs^(x u , x 2J . x 3J . x 4j ) co%^T^4cO 
?, tts^fi^- (y u , y 2 j. y 3 j)^. S54«ofT?iJo 

jfc^ajX§T'^fcff?iJcO#)fcfrm 1 , , m, 2> • ■ • 
£«oT:&#>4. 

[ 0 0 7 7 ] S 5 6<V^mi&TP I?£Xg{i. ^T<0 
-T4X*IT\ *&Tc9*S£{;L S5 3c9fi^tf(1«<0 

Tfc=5r4„ 

[0078] tLhcOUiH^-c*?^4 X ? lz . M^SIHI-C*) 



4ffig<0A*«#fc:^liE-r4^«l«-C&4aJ*fi^S:^ 
#>4i§^T\ #l)?l^-MiaF i;j coA5:fiEfflf 4^(i. 
«^#tt<^f«i-g;eW;:5&64 CI i: *rc# 4tf«\ 3£ 

4. 

[0079] ifcC *HBH<0feg^tt¥SI^^4 

rii, ( 3 ) 5^s;MW<7)ffiifew i;j t lt , mi <r>m 

^LA^fI^^*3(t4^-^ U y KERI^mfeP- 
HRFiji:, »2c7)^SScO»®-CfflV>^7 r -^^gM^ 
Ji ?:fflV^M?:St. 

[ o o 8 o ] z z x\ m-mm&A-ftmnt AtsiM^ 

(8) sSifcli (9) SCkRtTft"). ^-^SMSfe 
Ji ±^com2cOHtftcO^SST^L3t (16)^ 
- (22) setHtT'fc4. 

[0081] M^!ffiiffi-c-S)4ffisoai^fi^^je-r4 

^«MtTfc4A^m^**«)4^{i, «MM£HfflRP 

^mmii. ( 3 > &<x>m*itmtiMWii ttx a 

2) S«W ni ^fflV^ftbOt: (8) ^*fc{± (9) 
s^tC^L/SWuj ^fflv^41iyM±. H2^fc{iia3^ 

ffiT'S>4ffi-«coA^«-^^lC-r4^}i!Wi-C^4ai^« 
#**»4*&K-3vrci»!W*. 
[ 0 0 8 2 ] H5«, *^tOfefE»»tt^aK'&com 

^zttmi-zi'MmT'fo&iiijjm^zxfr&iir&iztift 

4iD^0-CT5-^-f7a-^-v-h. H6»i, I«lt<T-* 

t ZM®. Lte%&izii , ±^co^ 3 c^fife?)J&<Si: 
Pith ^4. «^(iS6 lc73Atii^llT-^^«SX^ 
t"£>4* { , c*U4Hlt:*iH-4S l lfOAaj^Hx-^ 

[ 0 0 8 3 ] iKC S 6 2 

6 3c07 f -^^SMISStlt^ffiiXg(i. 116 fc^SftT 
^4cOT. WTH6S:fflV^TittBy1t-4. 



(15) 



mm*? 10-262157 



[0084] S 7 1 <oatfMffl»£l«*C, m5<DS 6 

mT-fz&mmkLxmmL. s7 2jjrrcAaj* 

[ o o 8 5 1 mz s 7 2 <Dn?>-mf<7>i&$kco?kfem 
x\ ( 3 > rtra^neoa^fttt^oflwwij^flcaj 

ftfc LTte. ( 8 ) ^i^ti ( 9 ) ^^L^M^Mffl 

< . -mtrnzfeifr sit #t § s . 

[0086]»:S7 3<O^J<7>J&frgtbXirC\ 
<3) ^T'^L7t1W&^A^M^fix-?a^*T?'J£ 

*0<7)*S$-ffl^-r . ffWOlS^mjl, m^, • • • 
[00 87] <W;:s7 4<?)^«MSajlgT\ (6)jS 

*fctt massive, aCTfflfifr&TOHS*&fitf- 

HffSHttfAafl-f- (x,j. x 2 j, x 3 j)T'$>&0 
T. ttl*«^-(yij. y,j. y 3 J) s 7 3 ofi^ijco 
jftdrtrajXS-CSftfcf^KO*!^™,!. m,2, • ■ • 

[ 0 0 8 8 ] &?)S 7 5^H^Tfl^XS{±, ±X(0 

ajx*£<ciitf. 

[0089] S 7 6<7)x-^WSW^fi^a5XgJi. 
(16) (19) ^lC*L^-f-^»SWI8c^ 

xg«. -r<T«oAaj^T-^wcov^7-'-^»« 
*^T^-&tt. s 7 6<7>?~?mm®$m*n 

tHXgCcM^-C. ftcOAffl^flT-^ttlcovvCT-? 



tt«PM»a*tiJ*tt9£U *T<9*£-tt. E5H 
fcttS s 6 4 ^M^SHl^filXg^Jitr. 
[009 1 ] m5c0S6 4^a>S6 8tX'<r)fmXM 

{HJiaiat'cb^T. S6 4<7)X-m*hcr>T\ S6 4^ 

[0092] S6 4^S^«i!t<7)S:ttJXgT'ii. (3) 
n «<7)g*tf tfcD^ifcWi j £jra5^& . Z<7> 
*^M8*»A*e#*0*C, (8) sfcfcfctt 
( 9 ) *t*tfc«T«Mrc**A*fi#i: A7)ffi-^0 

»J3)c**^^S«tS^A7Jfi-^SS^i3(t&J--^ 'J 

(24) 3tT^jftLfc*><0*fflVvHtfJt< . -KtofcSfc 

[009 3] ULh. 05<7)S6 4fcW(i. T-?ffim 

6 l*^S68iT'^aOX^Srigl»^i:(cJ;0. 

±<0t&BJ§T';bA>-g> J; a . -r - ? *»£fWf«*£te "3 « 
^•C*->T4. JS^-SiJffi'C-^-S.ffiEcOA^fi-^tcWiiB-f 

[0094] «±, *5tw<?)mi&<7>4^(mmzwML 
tztf, z<nmzi>, m6iz7jkLtz7-?mmmwtkm2 

[0095] 

0 1li*7-«tt. 1 0 2— 1 0 4{±l&7UT--7' 
{±^7-S^A7J^afc LT^^-^df-A'^ 1 0 1 

miz-?\^xmw^&. 

[00 96] 07 (C^"T v-XxATii, 

1 0 l(i-^i.^<X7t^^-H«^^BX^TRGB^ 
ra«om#S:m*-ri>. ^7-X^-r-^l 0 ld^Jii^S 

1 <55ux-7';U 10 2-104 Icfc^TfipgtfffSrib 
ixS. 1&7CT— 7";H 02— 1 04COSWJ4. a^<± 
RSWcJtCT-rSRGBfe^na^r-^S:. iXW7hy 

y^X^XJf-^^gPl 0 5-C'L* a* b' feS^^T— 



(16) 
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=*A>^1 0 lc7)RGB(0#^>"9-^Sc7)jt^^iE-r 
*7-X#-v-M 0 llc^H^A^ J: 
££R' , G' , B' &m tffifcr&i, <fc 3(^^5:5^- 
dfcT'&S. lKjtf-/;H0 2-104t'R, G, 

y^SJ 1 0 5 T&g(C«c# L & V ^ft^feSfH L* ' a*' b 
*'K^LTaj^t-&. 10 

Tt^feS^L* a' b* COHr'-^t. VhU-y?X 
V^y^gSl 0 50tii*c7)L*' a*'b*'f^S<?Dr— 
*£-St3lM>. -IftKioT. #7-X3-^l 0 1 

T-?tfAjl2ix&Zblz%&. 
[ 0 0 9 7 ] 08 *SMH<Ofefiai»ttTiW5rffi*3« 

&ijkmm<7)^M<7)irmt lx, f-^MHt^ 

[0098] *T » S 8 1 C9AtfJ:*JSl7 J -:?*fPffiiX§ 
T. &S>a>ttf>fe3r«l5£L-C*33?fe£igL' a* b» £ 

0 UZWtX-Mhit. *7-^tio 1A> 

^>ai^$^&RGBfeisisi^fcftsft#^n-s!j-rs. i 
fi^ ct o , x^\m<^mf-9 1 Lx^m^mL* a* 

b* «0r— 9b. #7— -X^-^l 0 l«0{iS^JcORG 

«^<7)*7-^°-y^-fe{±. A:*j<ofe£fsii*K;:;t^ 
Wfy = - X,,)/L 0 ) 2 



akowKfcitsflpr*^, 3isi± ioomioo oat 
[0099] ij ^~mm(om^^t s . tt^-sf-K: 

coy 4 x^^mt^mumco^x^-^y^mrnxm 

ZZti^^XoltZXtt&bX^. £fz. 

-ffi(^g*SRGBM^f- 9 lz =7 yyMzfkh i. 
olzLlz*). Wmhox^-ryizftWiLXRGB&gm 

Wi.RGBm< l Z7>?'£.lzm&£olZ?2>k£\.\ 
0. RGBfeSra^x-^Srff-iBIL/if*. 7V*— 
[ 0 1 0 0 ] <Xlz. S 8 2<^fi*tt«iW»l£/>'^;< 

MVttw^jt-fzmfetz. mat*. (8) , 

(10), (12) 3CC«S-*-4a»ft»-W)IW*. 121 

Tfc***" (27) ~ (29) s^i:L-cia^ri». -r^ 

(l„ , a 0 , b 0 ) ^Khm^^mmtcomL. 
(r 0 . g 0 . b 0 ) **aj*m-9aa<oa»fb^ja». 

[0101] 
[&8] 

+ (0. - Oj)/ao) 2 
&i)/6o) 2 ) p + 1) (27) 



W 2 tf = ^((((mn^-L^ + Cm^-ay)) 2 

+(m 13 (6 i - 6,)) 2 )/(«o) 2 

+((m 2l (L,- - Lj)f + (m 22 (a, - a,)) 2 
+(m 23 (6 i - 6,)) 2 )/(Go) 2 

+((m 31 (Li - L^f + (m 32 (a, - ay)) 2 

+(m 33 (ft, - & 3 )) 2 )/(Bo) 2 ) p + 1) (28) 



[0102] m*Wfw<5x—fkii* zzx-ii (l 

0 . a 0 , b 0 ) , (R 0 , G 0 , B 0 ) . pc0388H 
-f6i<7>T&&. iRcoBVtX'ttfr&JiLdlZ. (L 0 , a 



(29) 

o , b 0 ) b (R 0 . G 0 , B 0 ) «i. ffi#±#<&& 



(17) 
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/a xm^ wu* x^- : jy^mm^m<~rtn£ 

[0 1 0 3] {ajxtf. ftS^'^^-^^fflSrM^UTA 

^mmz^ww&b, y-fxtism^ 

[0104] S«Jc73e?iJi: LTti. ttTOttfc LT L-* a* 
b* fe^/5f-i/3yf-^> Mx.fcf, a* =b 
• = 0T'L» #^L^o£;bSi:d&&^SIffi£*£«L 

Wi$£t&z\bhX1±&. 

[ o i o 5 ] ?»= , s 8 3 nwmmmffiinx'. &m 
cvm^wmcomfetz , mm-? 6 &mr/uzi vxj»b* 

;H0 2~104««[^t:^L, mz^hVv? 

<7>^(i, ljxjcT— x/n 02— i o4coimzmt> 

ZfzhblZ, WAXf—lV, tthib. a' = b* =0 

X'L* tfofrki oosx-mmzmkLx^&tefwm. 

1 0 0/2 5 6gjK^tTfcWi' 
raeti^V. -Vh^ X-?X*y7'& 10 5 CO 

mk£&#>&£iP>i<z\±, e&ffix-imiz&MZ'hz < l 

fc^fcv^dlWWrfetfL- a* b* ofg^x-*. « 
[0 10 6liXtC. S8 4CD^«!lffigajXg-C\ S83 

<m-x'bh<7>x~%mh. ztiizx*), m-mmx-h 

t>yvAx>r-ivb^* a*b* ffy&^s&T-fizftf 
$RG Bfe^tCfcit£^j|!|ffi^#£>fT.S. 



[0107] ifcfc:. S 8 5<Dfe*5S®&ifcJgffl7-'-?ai 
tHXST. S84i XnTMXfttzim.eM^mm.b =f> 

mm<r>?-?Mfrt>. &m®<r>m®£ik#>z>t:Mz<£m 

?$Hi#. £<^XxAco^JITlftoiEB£Ao7^ 

ttR»* fcV^ JUffiT'flfc-tirS *K WAX b L 
• a* b* cO^^T-'-^ORGBfe^^^JiHi, 
Oifalt. &M±10 0%Zc\z.Z,7 ! -?tf$&tlX^& 

™i#£> § . $> «ot- ? itQunMimmfem r- 

RGBM^fit^OA^ 1 0 OXcv&WtzXiX 

[0 1081KC s 8 6 <ofejt&«$ifc&£x*iT\ n 

XA-l 0 2 — 1 0 4^SSc£. X*r-)Wf-9 
irtti, l<k7cx-^>102~ 

io4i±, ummzitmt&RGB&gmn?-?*. 

BJJS5r#fi[tTSR' G' B' feSOOT'- 9 tz$mt 

stot-r*. r, g, B<7>&mm&Mznmt& i& 

TCx-T^l 0 2-1 0 4«0<Mfc<4, «Wil::R*;ttiG 
*>ttiB<^«Hltf>«£. «Wi{£L« c^ffiSri^T^V 

jat'^m-rsifc-RGBfe2gfa*^R' g' b* &s 

[0 109] -?h'; y^XV^^f^gpi0 5<0 

^SC^L* a* b' <0Mf-^^ffiot^6. 
L* a* b» <7)*S^7 ; -^(c:*tJE-r-S.RGBfeS^<0 
^SIHIS: . %(,zm>tz 1 <K7cr-X^ 1 0 2 - 1 0 4 £ 
Hit VCR" G' B' fe£ISkD-7*-*»C3SllU WR' 
G' B' fe^(7)f-^^VMJ 
0 5«A^]iL, A^3$ixS*9-ffi^L» a* b* 

&gkfficr)7-?*ntt<r>mmbLx. m^mk*®.^ 

^95 1 0 5 co«gc5r^-f & Cl t S . 
[0110] S 8 7<7)fe*QJi5i!2X;gT, -feja 

!SL)tV^fec7)L* a» b* fe^racOx-^tcML. feM 
Sf^c7)ai^T'^>-& LV a*' b* , feffi|B(?!)f r — 3* 

llfgL^ffiRcOfeL* a* b' Sra^ffifcLT^^ 
-X^f-^^ 1 0 l^ffi^T'^SRGBfe^P^x-^Sr 
^«!1L. -e#l$r 1 i^x-X/H 0 2-1 0 4£ffl^T 
R' G' B ' -feSraOx-^tC^IU Il:7f-ijy? 
XVX^fy^l 0 5tCj:-p-CL*' a*' b*'fe^^r 
-^fc^U. L* a* b« e&mztUf&T-?b(D 
&MiT$tt!>tl&£\,\ fflmcVZbZ. ffl«RGBf^ 

®<7)t-9 z&ft&b txmz-t&z b tMmx-fo h . 
zcobzu. &m<7)RGB&£m<7)T-?z& : f'mmb 



( 18) 
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• a' b* ^m^—^H-ML. RGB^S^f 
-?£l&7CX-7/H 02 — 1 04t-?hV y?X* 
VY^ 1 0 5 SrfflVvC L*' a*' b*'iZ&t&Lfzi><»t<7) 

[ 0 1 1 1 ] 09tt, *5feHB<^fe©H*tt^Hftfrat3i 
ffl-fS^rA<7)^2^?r^-tffi^;llT'i5&<, 0*. 
1 1 ltt3<jC7cr-yyK 112—11 5tel&7cx- 
?7k 1 1 6«i^5-7"Uy^-T'$>S. 
-M®tB^afcLTX?7-TUy^-l 16SrfflV\ 

^ximtz. 

[oii2j m9iz7pi-'>xf-j*T^ im-t^z*? 

#3frr— ://H 1 1 
1(2. L*' a*'b*'fe^g|c07-'-^^Y' M' C" K ' 

L* a* b* fe^ra^iittSr-^k^S^^-Str 
[0 113] 3&t£x-://H l l*^aJ^$i^SY , 

M' C" K' fe^ra^x-^f±. ^n^co^rtt 

1&7CT— 77H 12-11 5£AtJ$ix£. #1&7C 
tt-TOH 1 2-1 1 5l±. &f&ft<VMm&?f%^TY 
MCK&gffi&'f-fkLXXv-y'Vyf—l 1 6^ 
ftj]-t&« ZCOlttjCT—yjH 12—11 5 CO in w 
«2. 3^:76^-^1 1 lT'fi^T'^^rVUo^fe 

[0 114] *5-7'J^-l 1612. #l&7cx- 
:77H 1 2—1 1 5*>£>aj;fi£;ft6YMCKfe£|gJcDf 

* 7 -mm z mmi&mzzm-f h . 

c\<r>k£^ smf£tltz#y-WWt?>L* a' b' fe^fU 
tcfcft&ffifc.. 3»Scr-:77H 1 UcA^JSfLfc;*^ 
-■if- ^ <0 L * ' a * * b * ' M^Cfc ft 4 « fc * { -gC 

~$~h& 3 fe 3Wtf—^iV 1 l l coffin Jtf l &7cf 
— TVn l 2^1««i-|,ci:tc:<J;-?T. 5-x_t>1xtz*}? 

[ o n 5 ] H l ott. ^HB«efiBt«tt ; KIW6rife* 



xit, j>mn^kLx±M<nm2<7)9m<7)Bmxmw 

[ 0 1 1 6 ] lH\ S 9 1 OAffiflSr-^HlMKDS 
X\ YMCK-feS^COx— ^C0a^£0ffl^-€-C*^-/t 

YMCKfe^£0x-^tai*ffl9<0L* a* b* 

Aftflo&siap"] t t x^x-Sr < lx t > & ^ t a*a 

£L<. «;L»2\ YMCKW1 0%rhcO^T<7)ffl-^t 

0. YMCKfe^l^cOx-^^-fbtT^ajtJffllCOL* 
a « b* fi3aJ*^>7*-^*<**0SWfct=5rV^J:d=fe«« 
5:^|v^O^-4ii:(c:J;0. St 1 0 0*^ifcl 0 0 0 

< t>\,w#7-^-yj-$:iiij]-t&<o&mmmx'hz>. w 
iz, 1 6<n&imimtfm.mz^t 

Lx^hXo^mt. i. 'om^^mx-^y-^-/^- 

iz. tH^m^zit. ruv^-coHrtco^-tt^. a# 
mitt&'g.'&nffi.mz x h j * xtf^-> x \ »* tz*> . 

tr-a&L ±Mn*7-mmxftmm.<?>w&tmm&<F>X' 
[oii7] v&zs 9 2cr>m&tti-tmmmfi/^7*- 

*k*ww-^ftRj&r*. (9), (i 

1), (16)-(23) sSfcffiS-f 

«Tfc^ ( 3 0 ) - ( 3 9 ) ^cOJ; 3 Kt&rr 

tziz^mmoxwM^. <Yi , Mi , Ci , Ki > 

A^fi^OUT-*. (Y 0 , M 0 , C 0 , K 0 ) * 5 A 

^»m*jt{i^jaffi^ait»«^, , 3i , bi > 

WHttm^emr-f . (L 0 , a 0 . b 0 

[0 118] 
[»9] 



(19) #1^10-2 62 157 

+((Ci - C 3 )/C 0 ) 2 + ({A',- - KMKtfY + 1) (30) 

W 2ii = l/(((("»n(V;- - Yj)) 2 + (m 12 (Mi - Mj)) 2 

+(m 13 (C i - C,)) 2 + (mi 4 (K"i - Kj)) 2 )/(L 0 ) 2 
+((m 2l (Y i - IS)) 2 + (m 22 (Mi - M 5 )f 
+(m 23 (C,- - C s )f + (mu(Ki - A;)) 2 )/M 2 
+((m 31 (* - V,)) 2 + (ma(J4 - M,)) 2 
+(m 33 (C.- - C,)) 2 + (m 34 (A- i - A' i )) 2 )/(6o) 2 ) p + 1) 

(31) 

[$ti o] 



(20) 
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AL, = L'i — Li (32) 

Aai — a'i — at (33) 

Abi = b'i-bi (34) 

±Lu = (gflw*; - ^*») / (ijw*>) (35) 

A*. = (£<Vu*M - M)) / (£< V «*>) (37) 



W 3i = ((\(AL it + A£i) 2 x AL is x AL,| 1/4 /£ 0 
+|(Aa ia + Aa,) 2 x Aa j5 x Aai| 1/4 /a 0 
+|(A6 ia + A6i) 2 x Ab i4 x A&ip'VM" + 1) 
/((|(AL,-« - ALi) 2 x AL* x ALi\ l ^/L 0 
+|(Aa ia - Aa,) 2 x Aa is x Aa,| 1/4 /ao 
+|(A6.- a - A6,) 2 x Ab is x Afcf'W + 1) (38) 



[0 1 1 9] ZZr&X-ttWUJt-fMi. (Y 
0 , M 0 , C 0 , K 0 ) . ( L 0 , a 0 , b 0 ) . p<0 

o , M 0 , C 0 , K 0 ) t (L 0 , a 0 , b 0 ) Ji, ffl 

[ o 1 2 o ] <xfc. s 9 3cr>?-ftimmm*-wffiT. 

mX\ ( 1 6 ) - ( 2 3 ) j^cffiS-f i-%b*> 
( 3 2 ) - ( 3 9 ) 3«A»fea*W, 23i j s . * 



(39) 

[0121];XK. S 9 4 <?)^affi5*filiXST\ feJQ 

y/H 1 2-1 1 5c7)#£fc5r5fc£;tfe5£U <Xfc, 3<X7C 
^7^7^77-/^1 1 1^^^-I.t.^i:t 
I^tct— T'/H 1 2—1 1 5<omk 

zmb&fctbiz^ -r&jb^ m=c=k 

= 0%-C'Y*iO%*^ i o oxsTtymiz^kLT^z 

i>(T). Y=C=K = 0%T'M*«0%*>4.1 0 0%4T*% 
^tC^tUTV^tiO. Y = M=K = 0%-CC#0%*> 

hi ooxt-e&imiz^tix^zhc?). y=m=c = 



<2 1 ) 
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0%-CK**0%*>£> 1 0 0%iT-i^{C^LtT^S t 
** ; HWti:l/CRj&*-4. -fOXxyTfiDfllti. 
flaU{±'8t* y h^-y/U^ffio^i. 10 0/2 5 

[0 12 2] S&jcx-y/H 1 lcotMIScSr&tf) 

b* fc£OcOf-?fc. -ecO£K*tJfc-rSK?)%1I£«[ 

[0123] ifcfc:. S 9 5c^«!}ffigajI!IT\ S 9 4 

L» a* b' Q^CO?-? tRc^r-? |;j^t6Y 
MCfc£ficD^OIffi#f#£>:h.l> „ 
[ 0 1 2 4 ] <Xfc. S 9 6CDfe*^«f!rtK5£IgT\ SI 

;H 1 2-1 1 5<5DMfc£, jKO-r-^fctt 
-3T^-r*. Y, M. C, K<r>&$i.ft\iZfttf}T& 1<£ 

tcx-t'/h 12-11 5<r>wm.\3.. mmz&m&xy 
-fv^mmx-Jbh^ a* b* ^p a i<7)<a*^*aiL^ 

0%c7)feS^^-rSfelHt$:. «»KYifcttM4fcli 
C £MiKcO%m*b otf^^H^f-^ £T 
P7hL.7*D y br B K±!WilST'ifi(El^--S i fc tc J: 0 
Y" M ' C K' fe50A^YMCKfe?H^»W 

[0 12 5] ttlz, 3\^tr-^)Vl 1 10fi$lt£, 3 
ys^cx-y/H 1 10ttTjA£*fJ&? SL' a* b* 

t-tvh i inm-mznm-t . fjiti-j 

T3<y>fc Y M C fe£Ff3<7)%ffi i: T "ftC^T ifeSKO^ 
ffifcSr. 3fc£St«»t 1 <Xtct— 77H 1 2— 1 1 5 5rj£ 
WPKiiCAoTY' M' C K* feggHtf>-r-*fc: 

[0 126] fcrtfU £a#tti±3<fc7C?— 7*1 1 1 
tfOt&^^L*'a* , b* , fe^rBl^»T-^A < , Til 
>-?-l 1 6*<S^-Ci=6fetil*lc7)^tCffii.Sa*M 

«FiRJf, ltfSc-f-TVH 1 2-1 1 5 £0 
-10 0%crMXftfct ZOXliKK. ^HfS-fi6oT3 
iKjrr— 77H 1 lWf&^^L*'a*'b*'<OffiS^T 

<7>ift»W>£-f 6 i: , Y' M' C ' K ' fe^iaox-^^ 
-g|5#0%*}1i£;ttil 0 0%^^^. -t%:*> 
fittftS-frftfeY' M' C K* fe»f-^* { 



8u ilSG a m u t KMtmiix&^imt 

x. mim%WF&<oi&$it-*&z\ttfX'$&tf. mm 
ammth. 

[0127] m&ZS 9 7<D&mm&$tf^2J3xX\ 

mZLtz^&MMXt><F>L*' a*' b*'fe^|§Ii0T-^t 
*fU l6frt,fttlZtlZ>L' a* 

b* ggmof-frfk' 3^-5T^SA>$-p/<SC:i:£- 
L^T. SSSL^ffiS^Lx'a*' 
b * ' fe^H5<7)x- 9 Sr 3 iKTcr-- 7/H 1 1 T*S& LT 
Y' M' C ' K ' fe^^x-^?r^, -£<91I£1& 
TCX-TVP 1 1 2 — 1 1 5£fflWtYMCKfe£ia<Dx 
-^fclSa»L. YMCKfeS^x-^^S^SIMtL 

t*7-ryy?-i 1 6*^aj^j$ni>L* a* b* 

fe^HO-r-^^SIU AflSfifcL*' a*' b*'fe£ 

[ o 1 2 8 ] mz. *m\<?>m3<mi&cr>&mzmmi 
tz9mmtLx. m-mm.ifi\<X7VT-9<r>^\z-o^ 
x, mtmi, mszm^xwmt ±a<oj:3fc:. 

tJO. l<^7cr-^/H 02— 1 04t, 7h'j7?X 

vxdf-y^i o 5<7>®®&m>zm**tzmwi,tz. 

iftciO fe<S^14^-fb L , 1 ftjt^-y"^ 10 2- 

i o 4<r>m®<v&zi$m®mmzmtbWrfm&t:m%.L 
xwmtz . 

[0129] @8(i. T-?m%.ffl$l$:m^%\<^i%r£t<r> 

-oxwmtz. *?ssi<oxmimT-f*mmjM 

X\ $>e>fr tfrWlZZmfe LXb&rv-x. *r-ivn . y 
M L . A^ft^Hx- ^ k LT ^)Hg^ L' t XJim^ 

<nm7-?tLx<nRGBe&miz&w&7 : -?k<7)m 
hl-M<r>ftmzti$pthtf. ffl^iii o*^iti o< h 

^tmrnxS) ^xvx^f 105^ 
Wfrmwm^xA-isyrmmx'mzzn.^^* o iz 

[0130]^CS8 2<OS»ft(t00^K/>''«?^- 

r#)i><7)T\ (8) Steffi sf*a»ft»tH«t. (4 



(22) 
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0 ) s^J: o KKJWitfJ: k\ Zc\X\ ( L j ) AW 

?mm<r>xiim^. (Li) vxiim^cmT-'? . 
<l 0 > #A*e#2OT*»*HfoojSBL PjwwBBtt 

W^. = If ((((Li — Lj)/L Q 

[0 1 3 lll^ftltWt? (L 0 ) . p 
<7)2a^^^E*UTfe'3> XA-^>^^cOS 

Cl<7)^±, I<fc7&r-:7VH 0 2-104fcvh 
'J -y ?^7X# y ^gg 1 0 5 co«i-£&5i LT £>£ 

[ o 1 3 2 1 jMc. s 8 3comi>wm.mmiMTi±. i 

^76^-^ 10 2—104 coflgR&^&^<ofl9K 

BJJSL* #0*>4> 1 0 0£T'W;:3^^fcLT 
^^Sktf)^:|ft^8!^i6^:iS5g-r.S». tWf^TWW. 
fFUtf 8t*«/hcO^-^£te5*§£-l2:, 100/25 

[0 1 3 3] mz. S8 4c0^8!Hfi«a}XS-CJ±, m<0 

colli— -C*>S. ^c07^££T1»3£*n&-r£. 
h RG Bfe^ratC*3(tS^-SJffiA s #^fL5 . 

[ o i 3 4 ] . s s 5 <7>&wmmmm?- 9m 

ffiWT-^St*^. mmco&®L&iklib2>tzMZ < £W% 

^^{^(t-S^SS^i. Oifejff. ££(£1 0 0%&jg 

mtfQfrt, 1 0 0%eogEBfcAr>T^?.T-:?*r£Siai 

[0135] S 8 6*>6^fl«!tiftJglST\ H 

■W) 1 <}C7&r-y/uoflak*ifcJE^4 . 1 W£r-7)V 
1 0 2-1 0 4<D§HOB6Wi, iISaROTsteitW 

TWlOSTL* a* b' feSaB^-r-^fciHJWik 
£ £^8t»£*<ii < 2r£ J: o fc, R»*Tltt:V VM<3* 

&^<n*f-?\,zm&m-z\\LX'hh. ttz. m2<D 

Stmt. A**~>-<7)RGBJz>y-com!&<r>m^$:ffiEt 
Ztzib. X*-r1-l l z7'U-(WAjlZii.tzbZizR' G' 



im i] 

2 ) p + l) (40) 

R, G. B#*l:*fj£-$-£ l<fcjrr—:OH 02 
-104c^$t(2, mmzR£t:l l ±G$.M±Bcr> : fWm 
Oil*. cOffiSrfco-C^P-f X^-;K7) 

ittRGBfeSf^f.R' G' B' fe^P^iO^ra 

mz%.*>. zco&mm&ZM.^it-t&ztizmikfei- 
[oi36] mmz. s 8 7<r>&mmmmiM.x'& 

X^httm^hZk^Wfo^X^h . ftoT. ffBL 

-fORGBfe^^ai^T-'-^S-^aJL, -Hi.£l&7C 
f— 7VH 0 2—1 0 4£fflWCR' G' B' fe^^<7) 

5Srfflv^L* a* b» ^m^-^^mt. mm 

[0137] JjLb^i 3 tc. 1 &^T-77U<?)ffi8W>.* 

«sts*fflv^r i»:7c^f-'-^a^ff=5ri.s. 

[0138] Cf^erWif-^WMWfcffiv^firv* 

v\ 0 1 0 ^-fi: d =5rJ&a3iS^J: oT 1 <X7Ut— 7" 
fr<m%l<nmiibW. L£ft%:oZti>X'Z&. 
[0139] ^ri>. £ ZVt&JHB***rHm. &J) 

m^t'fmmb lx ( 4 0 > ^T'^L^s^w^fiss 

W lid 0»£ffli^FSB8Sfcff«rofcj&», ^(ctii^fi 

&*^mmt Lx=m9mzfT*dar&. mz_i*M9£ 
7Frr*>?-Tv yf-o&o&m&tz&^xb, fr^o 

[0140] J: o C 
[0141] 



(23) 



#^¥1 0-2 6 2 1 57 



•o xm^mw * 7 ^ - ? * & ittios^f 
[Hi] *ira<7>feffie»« ; f«*i6 w»t i <nmM<7) 



[04 3 *mtofe&amfrm&&im 3 oiufe^ 
jB»te*iviT. mi-Mmx'Z>&fzmcoxj]m-tTiztt®-? 
^mmx-hhtatim^s^ ®&iz&v siwo- 

[05 3 *«BJcOfee^#tti^S!J^^m4<7)|ISt<0 

&?mmx'h h £j)m^z**>z ig&iztmmttn- 
[06 3 *ira<ofefiaa«^»srte^4<osaiw 

[073 *mi<tt&mm?mm*mm-t' h ^x 
y-^m 1 comz^-tffifmxb h . 
[08 3 *«w«feffii»tt^aft«£*affl-t*^ 
vs^se 1 & fesaa^i^^fe x x/mssm 

CO— #J£jjrf 7 C7— f-v— h T'3fe h . 

[093 *IMBOfeej*l»tt^«t)Sr86*3WW-* J'X 
•r Acom 2 <O0iJ£Srf fl§j£0T& & . 

[0103 4^<0fifi3M^5 ; attfi£*affl-r* ^ 

XT-^com 2 oeajtctj itsfe^a^ifc^fc itMiig® 
[**f-«ift(y§3 

10 1-*7-^tt. 10 2~104---l<5:7C7 i - 
1 0 5 ■■••7hU ytX^X^yftt. 111—3 
ftjCT—7fo, 1 1 2~1 1 5- • l&TC^-TVK 1 1 



[07 3 



[01 03 




a 

b* 





G 




G' J 





















S 9 1 

1 xmiimr-fMrn^jim 1 



S92 



1 



S93 



S 94 



I 



S95 



I 



S96 



1 S9 7 
-J^ 

<3DD 



(24) 



0-2 62 157 



[HI] 

S 1 1 

I AaiA%f-^f«ii8 1 
1" S 1 3 

j S15 

I ?«siigax# 

-j sTtT 

j sn 




I 



(25) 



#gfl¥l 0-2 6 2 1 5 7 



[02] 



S 12 



S 13 



S14 



S15 



S16 



s 17 



S 18 



S 1 9 




SI 1 

1 Atiit>nT-mmmi.n \ 




J 



(26) 



ftm^l 0-2 6 2 1 57 



IH33 



S3 1 



S32 
1 S3T 

1 sI7 





(27) 



10-262157 



[H4] 



[05] 



J~ S5 1 



S52 



S53 



1 S54 , 





S 6 1 




**T3MifSx*s 





1 S66 

1 ff?iJ(Plfc»j?HJXjg 



7 S6 7 

1 ?iiimxs i 



*!T¥J£xfg 



S68 



NO 



(28) 



10-262157 



CH6] 



[H8] 




S 7 1 



S72 



i 



S73 



S74 




S 7 6 



S7 7 



S81 



|xm^giT-^M^iii§xe| 

[ S82 



I 



S83 



1 S 8 4 

1 S 8 5 



ataxia 



i 



S 8 6 



1 &mm&&.wmx.& 



I 



Yes 



(29) 



%fffl¥ 10-262 1 57 




Searching PAJ 



Page 1 of 2 



PATENT ABSTRACTS OF JAPAN 

(1 1 )Publication number : 1 0-2621 57 

(43)Date of publication of application : 29.09.1998 



(51)lnt.CI. 




H04N 1/60 
H04N 1/46 


(21 Application number 


09-332717 


(71 Applicant : FUJI XEROX CO LTD 


(22)Date of filing : 


03.12.1997 


(72)lnventor : IKEGAMI HIROAKI 



(30)Priority 

Priority number : 09 4949 Priority date : 14.01 .1997 Priority country : JP 



Sit 



(54) METHOD FOR PREDICTING COLOR TRANSFER CHARACTERISTIC 
(57)Abstract: 
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signal and to predict an input signal with respect to an 
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input signal from part of the input signal and an optional 
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SOLUTION: Input/output real data pair are prepared 
(S1 1) and a predicted value is prepared in the S12. An 
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using the data, and a component of a matrix and a 
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CLAIMS 



[Claim(s)] 

[Claim 1] It connects with a matrix so that it may become the linear relation characterized by providing 
the following. A component of this matrix is determined that the square sum of a value which carried 
out weighting to difference for every signal component of an output forecast calculated using this matrix 
and two or more output live data corresponding to it by coefficient of weighting of forecast-ed 
dependence will become min from live data of two or more input signals. A color transfer-characteristics 
prediction method characterized by calculating an output forecast using this matrix from an input signal 
which is a forecast-ed Live data of two or more input signals of a color picture input unit or a color 
picture output unit It sets to a color transfer-characteristics prediction method of searching for an output 
signal which is a forecast corresponding to an input signal of arbitration which is a forecast-ed from a 
data pair of live data of an output signal corresponding to it, and is a constant term about relation 
between an input signal and an output signal. 

[Claim 2] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast-ed, and an input signal by constant for 
every input signal component — difference — a color transfer-characteristics prediction method 
according to claim 1 characterized by being the monotonically decreasing function of Euclidean distance 
in standardization input signal space for which it asks from a component. 

[Claim 3] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast-ed, and an input signal by constant for 
every input signal component — difference — with a monotonically decreasing function of Euclidean 
distance in standardization input signal space for which it asks from a component A color transfer- 
characteristics prediction method according to claim 1 characterized by consisting of data precision 
functions calculated from an interrelation of an output forecast beforehand predicted to live data of each 
input signal, and live data of an output signal. 

[Claim 4] It connects with a matrix so that it may become the linear relation characterized by providing 
the following. So that the square sum of a value which carried out weighting to difference for every 
signal component of an output forecast calculated using this matrix and two or more output live data 
corresponding to it by coefficient of weighting of forecast-ed dependence and component dependence of 
this matrix may become min from live data of two or more input signals A color transfer-characteristics 
prediction method characterized by determining a component of this matrix, a coefficient of weighting, 
and a forecast by the successive approximation method Live data of two or more input signals of a color 
picture input unit or a color picture output unit It sets to a color transfer-characteristics prediction 
method of searching for an output signal which is a forecast corresponding to an input signal of 
arbitration which is a forecast-ed from a data pair of live data of an output signal corresponding to it, and 
is a constant term about relation between an input signal and an output signal. 

[Claim 5] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast-ed, and an input signal by constant for 
every input signal component — difference — with the 1st monotonically decreasing function of 
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Euclidean distance in standardization input signal space for which it asks from a component By dividing 
by constant for every signal component, after changing into a component of output space which 
considered sensitivity using a component of this matrix by difference for every signal component with 
live data of an input signal of a forecast-ed, and an input signal A color transfer-characteristics 
prediction method according to claim 4 characterized by consisting of the 2nd monotonically decreasing 
function of Euclidean distance in standardization output signal space for which it standardizes and asks. 
[Claim 6] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast-ed, and an input signal by constant for 
every input signal component — difference — with the 1st monotonically decreasing function of 
Euclidean distance in standardization input signal space for which it asks from a component The 2nd 
monotonically decreasing function of Euclidean distance in standardization output signal space for 
which it standardizes and asks by dividing by constant for every signal component after changing into a 
component of output space which considered sensitivity using a component of this matrix by difference 
for every signal component with live data of an input signal of a forecast-ed, and an input signal, A color 
transfer-characteristics prediction method according to claim 4 characterized by consisting of data 
precision functions calculated from an interrelation of an output forecast beforehand predicted to live 
data of each input signal, and live data of an output signal. 

[Claim 7] The color transfer-characteristics prediction method characterized by to determine the 
component of this matrix, the coefficient of weighting, and a forecast by the successive approximation 
method so that the square sum of the value which is characterized by to provide the following, and 
which carried out weighting by coefficient of weighting of forecast dependence at least to difference for 
every signal component of an output forecast which tied up with a matrix so that it might become linear 
relation, and was calculated using this matrix from live data of two or more input signals, and two or 
more output live data corresponding to it may become min Live data of two or more input signals of a 
color picture input unit or a color picture output unit It sets to a color transfer-characteristics prediction 
method of searching for a part of input signals which are a forecast corresponding to an output signal of 
arbitration which is a forecast-ed, or remaining input signals which are the forecasts corresponding to a 
part of an output signal and an input signal of arbitration which is a forecast-ed from a data pair of live 
data of an output signal corresponding to it, and is a constant term about relation between an input signal 
and an output signal. 

[Claim 8] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast, and an input signal by constant for 
every input signal component - difference - a color transfer-characteristics prediction method 
according to claim 7 characterized by being the monotonically decreasing function of Euclidean distance 
in standardization input signal space for which it asks from a component. 

[Claim 9] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast, and an input signal by constant for 
every input signal component - difference - with a monotonically decreasing function of Euclidean 
distance in standardization input signal space for which it asks from a component A color transfer- 
characteristics prediction method according to claim 7 characterized by consisting of data precision 
functions calculated from an interrelation of an output forecast beforehand predicted to live data of each 
input signal, and live data of an output signal. 

[Claim 10] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast, and an input signal by constant for 
every input signal component - difference - with the 1st monotonically decreasing function of 
Euclidean distance in standardization input signal space for which it asks from a component By dividing 
by constant for every signal component, after changing into a component of output space which 
considered sensitivity using a component of this matrix by difference for every signal component with 
live data of an input signal of a forecast, and an input signal A color transfer-characteristics prediction 
method according to claim 7 characterized by consisting of the 2nd monotonically decreasing function 
of Euclidean distance in standardization output signal space for which it standardizes and asks. 
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[Claim 11] input space standardized when a coefficient of said weighting divided a difference for every 
signal component with live data of an input signal of a forecast, and an input signal by constant for 
every input signal component — difference — with the 1st monotonically decreasing function of 
Euclidean distance in standardization input signal space for which it asks from a component The 2nd 
monotonically decreasing function of Euclidean distance in standardization output signal space for 
which it standardizes and asks by dividing by constant for every signal component after changing into a 
component of output space which considered sensitivity using a component of this matrix by difference 
for every signal component with live data of an input signal of a forecast, and an input signal, A color 
transfer-characteristics prediction method according to claim 7 characterized by consisting of data 
precision functions calculated from an interrelation of an output forecast beforehand predicted to live 
data of each input signal, and live data of an output signal. 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[The technical field to which invention belongs] This invention reads a color copy, performs an image 
processing, and relates to the color transfer-characteristics prediction method of color picture I/O 
devices, such as the scanner and printer which are used in case the color processing coefficient of the 
various image processing systems performed in order to obtain the color reproduction which wishes to 
have a manuscript image in the digital full colour copying machine reproduced on a record medium-ed, 
color facsimile, a color system, etc. is determined, and a display. 
[0002] 

[Description of the Prior Art] For example, after performing the image processing which changes a 
manuscript into the signal of the color space for which it does not depend on equipment from the signal 
of reading and the RGB color space which is the output of a scanner with a scanner and carrying out a 
certain edit processing, the image processing changed into the signal of the CMYK color space which is 
the input of a printer is performed, and a system which carries out a printed output by the printer is 
assumed. In this case, from the signal of the RGB color space depending on the equipment which is the 
output of a scanner, performing the image processing changed into the signal of the color space 
independent of equipment will grasp correctly the relation between the color of an input of a scanner, 
and the color of the RGB color space of the output of a scanner, i.e., the color transfer characteristics of 
an input unit, and it will perform the image processing of that reverse. Performing the image processing 
which similarly is changed into the signal of the CMYK color space depending on the equipment which 
is the input of a printer from the signal of the color space independent of equipment will grasp correctly 
the relation between the color in the CMYK color space of an input of a printer, and the color of the 
output of a printer, i.e., the color transfer characteristics of an output unit, and it will perform the image 
processing of the reverse. Even if an output unit is a display, the same thing can say. 
[0003] It is a base to grasp the property of equipment correctly, although the methods of deciding the 
color processing coefficient used in an actual image processing differ by whether how a difference of an 
image-processing method and the color reproduction to wish set up, and if predict the output signal over 
the input signal of the arbitration of a color-picture input unit or there is a method which has predicted 
the input signal over the output signal of arbitration enough, and can do it, it will become easy to decide 
the color processing coefficient of an image processing. 

[0004] There are what predicts the physical property of an I/O device as the color transfer-characteristics 
prediction method of a color I/O device using the physical model grasped and modeled, and the method 
of processing the live data of I/O of equipment statistically and predicting them. 
[0005] For example, in JP,5-18305,B, it predicted by having assumed the physical model called a noy 
violence UA equation in 3x3 matrices in the printer which is an output unit in the color scanner which is 
an input unit, and the technique of asking for a masking parameter by convergence count is proposed. 
However, even if, as for such a physical model, model precision determines the color processing 
coefficient of an actual I/O device well, what can be equal to practical use is not obtained. Moreover, 
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development of a physical model with more high predictability had the problem that it was necessary to 
devise a model for every equipment of the, when the classes of equipment differed, even if it is difficult 
and development of a model was completed. 

[0006] On the other hand, between I/O of a printer, the method currently generally called computer color 
matching assumes a high order polynomial or the function corresponding to the physical property of 
equipment, determines the coefficient of a function from the data pair of the live data of two or more 
input signals, and the live data of the output signal corresponding to it, and predicts an input signal [ as 
opposed to the output signal of arbitration for the function ] by solution Lycium chinense 
asymptotically. However, predictability was influenced from the relation of the configuration of the 
function for which live data took and it depended on equipment the direction, and when it was the high 
order polynomial which requires computation time, there was a problem of being able to extrapolate 
outside of the color space of live data and being unable to predict it. 

[0007] Moreover, from the data pair of the live data of two or more input signals, and the live data of the 
output signal corresponding to it, using the neural network, to P. 125-129, the output signal over the 
input signal of arbitration is predicted, or the method of predicting the input signal over the output signal 
of arbitration is proposed in "high-degree-of-accuracy printer model by the high-degree-of-accuracy 
color conversion-neural network by flexible GCR -" besides JP,7-87347,A and Murai, the Society of 
Electrophotography of Japan, and Vol.35 No. 2 or 1996 years. By this method, when study took time 
amount too much or the outside of the color space of live data was extrapolated and predicted as well as 
use of a polynomial, there was a problem that it could not be used. 

[0008] Moreover, in JP,2-226870,A, the method of tying up 3 local **** space of I/O with a linearity 
matrix, and predicting the input signal over the output signal of arbitration, when the I/O signal after 
carrying out, the increase of the number of data pairs and is a three dimension is proposed by 
interpolation from the data pair of the live data of two or more input signals, and the live data of the 
output signal corresponding to it. By this method, since the smoothing feature of data is not contained, 
when the live data of an I/O signal contain a noise, it cannot be coped with. Moreover, this method is the 
technique which can be used when the dimension of an input and an output is the same dimension, and, 
in the case of the three dimension as which an output expresses a color in four dimensions [ like 
YMCK ] whose input is, cannot cope with it. Furthermore, when the outside of the color space of live 
data was extrapolated and predicted, there was also a problem that it could not be used. 
[0009] Moreover, on the U.S. Pat. No. 5471324 specifications, the method of deciding the coefficient of 
the conversion look-up table for printers is proposed by simple interpolation by predicting an input 
signal from an output signal with the weighting average the increase of the number of data pairs, and 
after carrying out from the data pair of the live data of two or more input signals of a printer, and the live 
data of the output signal corresponding to it. However, this method also had the almost same problem as 
above-mentioned JP,2-226870,A. 

[0010] Moreover, the method of determining the coefficient of color processing as JP,2-289367,A using 
the technique guessed are statistical from the data pair of the live data of two or more input signals and 
the live data of the output signal corresponding to it is proposed. However, it is expected that there is a 
problem that the continuity of the forecast in the boundary of subspace is not secured though description 
of interpolating and concreteness being missing only by being described by the language of**, and 
there being the same problem as above-mentioned JP,2-226870,A, or dividing into subspace to the 
smoothing feature which divides the contents into subspace and which is made to correlate is contained. 
[0011] 

[Problem(s) to be Solved by the Invention] How to have made this invention in view of the situation 
mentioned above, process the live data of I/O of a color I/O device statistically, and predict color 
transfer characteristics, Namely, the method of predicting the output signal over the input signal of 
arbitration and the method of predicting the input signal over the output signal of arbitration, Or it aims 
at offering a method with the following features in a concrete form as a method of predicting a part of 
remaining input signals from a part of the output signal and input signal of arbitration. 
1 . It is a method independent of the class (physical characteristic) of equipment. 
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2. It is the method which has a smoothing feature on the assumption that the live data of I/O contain the 
noise. 

3. Even when the dimension of an input signal is four dimensions, treat on a par with the case of a three 
dimension. 

4. Predictability is high. 

5. The continuity of a forecast should be secured. 

6. There needs to be extrapolation capacity. 

7. As compared with the method using the conventional neural network, the time amount of prediction is 
short. 

[0012] 

[Means for Solving the Problem] This invention ties up relation of an input signal and an output signal 
of a color picture input unit or a color picture output unit with a matrix so that it may become linear 
relation including a constant term, and it determines a component of this matrix from a data pair of live 
data of two or more input signals, and live data of an output signal corresponding to it. At this time, an 
input signal is made into a forecast-ed, an output signal is made into a forecast, and a component of a 
matrix can be determined if it is made for the square sum of an output forecast calculated using a matrix 
from live data of an input signal and a value which carried out weighting to difference for every signal 
component with two or more output live data corresponding to it by coefficient of weighting of forecast- 
ed dependence to become min. Moreover, what is necessary is just to determine a component of a 
matrix, a coefficient of weighting, and a forecast by the successive approximation method so that the 
square sum of a value which carried out weighting to difference for every signal component of an output 
forecast and output live data by coefficient of weighting may become min in using a coefficient of 
component dependence of a matrix with forecast-ed dependence as a coefficient of weighting. 
[0013] An output signal is made into a forecast-ed and an input signal is made into a forecast. An output 
signal and some input signals Or a forecast-ed, So that the square sum of a value which made the 
remaining input signals a forecast and carried out weighting to difference for every signal component of 
an output forecast calculated using this matrix from live data of an input signal and two or more output 
live data corresponding to it by coefficient of weighting of forecast dependence may become min A 
component of a matrix, a coefficient of weighting, and a forecast can be determined by the successive 
approximation method. In this case, a coefficient of component dependence of a matrix may be used 
with forecast-ed dependence as a coefficient of weighting. 

[0014] Moreover, in any case, precision may be raised, using as a component a data precision function 
calculated from an interrelation of an output forecast beforehand predicted to live data of each input 
signal as a coefficient of weighting, and live data of an output signal. 
[0015] 

[Embodiment of the Invention] About the gestalt of operation of the 1st of the color transfer- 
characteristics prediction method of this invention, general explanation is performed first. Although an 
input signal explains here using (7) types from (1) type which takes for an example the case where an 
output signal is a three dimension in a three dimension or four dimensions, and is shown below, the 
number of dimension of an I/O signal is not necessarily limited to these, and is the same at other number 
of dimensions. (1) By the three dimension, as for a formula and (6) types, the input signal corresponds, 
when an output signal is a three dimension, (2) types and (7) types correspond, when an input signal is 
four dimensions, and in both cases, (3) types correspond. 
[0016] 
[Equation 1] 
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[0017] (1) type and (2) types First, the live data (xli, x2i, x3i) of the input signal of an n-tuple, or (xli, 
x2i, x3i, x4i) i= 1 - n, connecting with a matrix is shown so that it may become the relation of linearity 
including a constant term about the forecast (yli\ y2i f , and y3i - i= 1 - n, and it has become each 
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component of a matrix of ml 1, ml2, and ... 

[0018] The following (3) types are the square sum Ej of the Euclidean distance of the output forecast 
calculated using the matrix shown in (1) type or (2) types from the live data of two or more input 
signals, and two or more output live data corresponding to the live data of the input signal by which 
weighting was carried out. It is shown, here — (-- y — one — i — ' — y — two — i — 1 — y ~ three — i — 1 — ) 

- i __ = „ one - n — (1) — the live data (xli, x2i, x3i) of the input signal of an n-tuple shown by the 
formula or (2) formulas — or (xli, x2i, x3i, x4i) It is the output forecast calculated using the matrix from 
i=l-n. (yli, y2i, y3i), and i= 1 - n It is live data of the output signal of the n-tuple corresponding to i= 1 - 
n. (xli, x2i, x3i) Or (xli, x2i, x3i, x4i) Wij, i= 1 - n It is the coefficient of weighting to the output 
forecast calculated using the matrix shown in (1) type or (2) types from the live data of these input 
signals, and the Euclidean distance between the live data of an output signal. 

[0019] (6) As for the formula, how an input signal calculates a forecast from the forecast-ed of 
arbitration in case an output signal is a three dimension in a three dimension is shown, and ml 1, ml2, 
and ... are each component of the same matrix as (1) type. When a forecast-ed is an input signal (xlj, 
x2j, x3j), the output signal (ylj, y2j, y3j) which is a forecast easily can be searched for by substituting 
for (6) types. When a forecast-ed is an output signal (ylj, y2j, y3j), the input signal (xlj, x2j, x3j) which 
is a forecast by solution Lycium chinense about (6) types at reverse can be searched for similarly. 
[0020] (7) How, as for a formula, an input signal calculates a forecast in case an output signal is a three 
dimension in four dimensions is shown, and ml 1, ml 2, and ... are each component of the same matrix as 
(2) types. When a forecast-ed is an input signal (xlj, x2j, x3j, x4j), the output signal (ylj, y2j, y3j) 
which is a forecast easily can be searched for by substituting for (7) types, the input signal which is an 
input signal which it specifies as the part, for example, (x4j), the forecast-ed, of an input signal, and (7) 
types are solved conversely, and is the remaining forecast since (7) types are not made as for solution 
Lycium chinense to reverse when a forecast-ed is an output signal (ylj, y2j, y3j) — for example, (xlj, 
x2j, x3j), it can ask. 

[0021] Thus, even when the dimension of an input signal is four dimensions, it can treat on a par with 
the case of a three dimension. Of course, irrespective of a number of dimension, it can predict and the 
application range is wide. Moreover, since input / output relation is tied up with comparatively simple 
linear relation, extrapolation capacity is strong, even if the live-data pair of I/O does not cover the whole 
region of an I/O color space, in the whole region of an I/O color space, it is usable, and it can be used 
also for the property prediction besides the color gamut of an I/O color space. 
[0022] By using the well-known method called the least square method, when the coefficient Wij of 
weighting, i= 1 - n have become settled, it is Ej. It can ask for each components mil and ml 2 of a 
matrix, and ... on conditions which are made into min. However, since it is not uniquely decided when 
dependent on the input value the coefficient Wij of the aforementioned weighting, i= 1 - whose n are the 
components or forecasts of a matrix, it is Ej of (3) types. Under the conditions of considering as min, the 
coefficient Wij of weighting, the components ml 1 and ml2 of i= 1 - n, and a matrix, and the optimum 
value of ... and a forecast are determined using the technique of successive approximations. 
[0023] In addition, Ej Since it can decompose into the sum of a positive component about each signal 
component of output signal space as shown in (4) types which transformed (3) types, considering as min 
is equivalent to it being independent and making into min the amount decomposed into each signal 
component, the difference of two or more output live data corresponding to the output forecast and it 
which asked for the following (5) types using the matrix from the live data of two or more input signals 

— every signal component - constant twice — square sum Ej' of the distance by which weighting was 
carried out to the value carried out is shown. The method of this invention can completely be similarly 
applied, when making Ej 1 as shown in (5) types into min. 

[0024] moreover, the input space where the coefficient Wij of weighting standardized the difference for 
every signal component with the live data of the input signal of a forecast-ed or a forecast, and an input 
signal — difference - with the monotonically decreasing function Fij of the Euclidean distance in the 
standardization input signal space for which it asks from a component the monotonically decreasing 
function Gij of the Euclidean distance in the standardization output signal space for which it 
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standardizes and asks after changing into the component of the output space which considered 
sensitivity using the component of this matrix by the difference for every signal component with the live 
data of the input signal of a forecast-ed or a forecast, and an input signal ~ since — it constitutes. Thus, 
by constituting the coefficient Wij of weighting by the monotonically decreasing function of Euclidean 
distance, about live data with the large difference of distance, weight can be made small, effect can be 
lessened, weight can be enlarged about live data with the small difference of distance, and it can treat as 
important data. Moreover, since these function configurations are monotonically decreasing functions, 
the continuity of a forecast is secured theoretically, and the coefficient of color processing of a color 
picture I/O device can be decided without caring about the discontinuity between parts. 
[0025] An output signal takes for an example the case where an input signal is a three dimension in a 
three dimension or four dimensions about this, and explains using (12) types from (8) types shown 
below. (8) A formula and (10) types correspond, when an input signal is a three dimension, when an 
input signal is [ the output signal of (9) types and (11) types ] a three dimension in four dimensions, 
correspond, and in both cases, (12) types correspond. 
[0026] 
[Equation 4] 
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+((*3i - X 3 j)/x 3 o) 2 ) 
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Wuj = F {j (((x li -x xj )ix lQ y + ((x 2i -x 2i )ix 2{i y 

+{{x*i - X 3 j)/x Z0 ) 2 + ((^4i - Xij)/X i0 ) 2 ) 
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W 2ij = Gy(((m u (»u - X!,)) 2 + {m 12 (x 2i - x 2j )) 2 
+(m 13 (a;3 { - x 3j )f) j {y lQ f 
+({m 2l (x u - xij)) 2 + (m 22 (x 2i - x 2 j)) 2 
+(rn 2Z {x zi - x 3j )) 2 ) / (y 20 f 
+((m 3 i(x li - a? y )) 2 + (m 32 (x 2i - x 2j )f 



W 2ij = Gijdimnixu - x X j)) 2 + {m l2 {x 2i - x 2j )f 

+(m 13 (a; 3i - x zi )f + (m 14 (x 4i - x 4j )) 2 )/ {y 10 f 
+((m 2i (x u - xxj)) 2 + {m 22 (x 2i - x 2j )) 2 

+(m 23 (^ 3 i - X3j)f + (m 24 (x4i - ^4;)) 2 )/(y20) 2 
+(( m 3l( X U - X lj)Y + ( m 32(x 2i - X 2 j)f 

+(^33(2;^ - x zj )f + (m 34 (x 4i - x 4j )) 2 )/ (y 30 ) 2 ) (11) 
W l2ij = H(W lih W 2ij ) (12) 



[0027] first, the input space which standardized the difference for every signal component of the input 
signal of the forecast-ed or forecast (8) types and whose (9) types are the elements of the coefficient of 
weighting, and live data — difference — the monotonically decreasing function of the Euclidean distance 
in the standardization input signal space for which it asks from a component is shown. Here (xlj, x2j, 
x3j) (xlj, x2j, x3j, x4j) Or the input signal of a forecast-ed or a forecast, (xli, x2i, x3i) Or (xli, x2i, x3i, 
x4i) the live data of an input signal or (xlO, x20, x30) (xlO, x20, x30, x40) the constant of 
standardization of input signal space, and Fij are monotonically decreasing functions. By (8) formulas, 
for example, the difference for every signal component with the live data (xli, x2i, x3i) of the input 
signal (xlj, x2j, x3j) of a forecast-ed or a forecast, and an input signal (xli-xlj), It asks for (x2i-x2j) and 
(x3i-x3j), standardizes by the constant (xlO, x20, x30) of standardization of this, and carries out to (xli- 
xlj) / xlO, (x2i-x2j) / x20, and (x3i-x3j) / x30. These square sums are made into the parameter of a 
monotonically decreasing function Fij. 

[0028] Next, (10) types and (1 1) types show the monotonically decreasing function of the Euclidean 
distance in the standardization output signal space for which it standardizes and asks, after changing the 
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difference for every signal component with the live data of the input signal of a forecast-ed or a forecast, 
and an input signal into the component of the output space which considered sensitivity using the 
component of this matrix. The input signal of here, a forecast-ed, or a forecast, (xli, x2i, x3i) Or (xli, 
x2i, x3i, x4i) the component of the matrix of the above-mentioned [ the live data of an input signal ml 1 
and ml2, and ... ], the constant of standardization of (ylO, y20, y30) of output signal space, and Gij are 
monotonically decreasing functions. After hanging the component of the difference for every signal 
component with the live data of the input signal of a forecast-ed or a forecast, and an input signal and 
matrix by changing into the component of the output space which considered sensitivity here and 
carrying out a square for every signal component, compared with the method of meaning taking the sum, 
taking the sum, without carrying out a square, and changing into the Euclidean distance of the usual 
output space, the difference is large. When taking the sum, without carrying out a square, even if an 
absolute value is large, depending on the sign of each item, the sum may become small. (10) The 
absolute value of each term can be considered as sensitivity by carrying out a square like a formula and 
(11) types. 

[0029] For example, (10) types ask for the difference (xli-xlj) for every signal component with the live 
data (xli, x2i, x3i) of the input signal (xlj, x2j, x3j) of a forecast-ed or a forecast, and an input signal, 
(x2i-x2j), and (x3i-x3j). Although it is convertible for the component of output space by imposing a 
matrix on this, the sum is taken after hanging and carrying out the square of the component of a matrix 
as mentioned above here, that is () [ mil] (xli-xlj) 2+(ml2 (x2i-x2j)) 2 + 2, 2 (m21 (xli-xlj))+(m22 
(x2i-x2j)) 2+(m23 (x3i-x3j)) 2, and 2 (m31 (xli-xlj))+(m32 (x2i-x2j)) 2 (ml 3 (x3i-x3j)) + (m33 (x3i- 
x3j)) 2 It asks. Using the constant (ylO, y20, y30) of standardization, the division of these is done, they 
are standardized, are added by 2, 2 (x20), and 2 (x30), respectively (ylO), and it is considering as the 
parameter of a monotonically decreasing function Gij . 

[0030] Next, it is shown that (12) types compound two above-mentioned monotonically decreasing 
functions, and create the coefficient of weighting. Here, it is necessary to compound the method of 
compounding for example, so that the feature that two functions which are the element are 
monotonically decreasing functions of each Euclidean distance may not be broken down, and an 
operation like the sum or a product can be used for it. 

[0031] Thus, W12ij calculated by (12) formulas is the coefficient Wij of weighting of (3) types. Ej 
shown in (3) types with the least square method using the coefficient Wij of this weighting It asks for 
the components ml 1 and ml2 of the matrix made into min, and ... And a forecast is calculated using the 
matrix which consists of components ml 1 and ml 2 of the matrix searched for, and ... Furthermore, if 
needed, the coefficient of weighting is again calculated using this calculated forecast, and it is Ej. It asks 
for the components ml 1 and ml2 of the matrix made into min, and and a forecast is re-calculated 
using this. By repeating such processing, it can converge and a forecast can obtain a desired forecast. 
[0032] Drawing 1 is a flow chart which shows an example of the actuation in the gestalt of operation of 
the 1st of the color transfer-characteristics prediction method of this invention. The live data (xli, x2i, 
x3i) of the input signal of the n-tuple which first was stated by (1) - (3) formula [ two or more sets of 
live-data pairs of an I/O device which want to actually predict a property at the I/O live-data pair 
preparation production process of Sll, i.e., the above-mentioned, ], and i= 1 - n, The live data (yli, y2i, 
y3i) of the output signal of the n-tuple corresponding to it, i= 1 - n are prepared. As the example, in the 
case of a color scanner, the various color patches which have measured the color beforehand are made to 
read, the RGB value of an output is measured, and it is the chromaticity coordinate of an input, for 
example, L* a* b*. It is equivalent to preparing a live-data pair with the RGB value of an output. 
Moreover, in a display, it is the RGB value of an input and the color coordinate of an output, for 
example, L* a* b*, by displaying a color patch in the various combination of the RGB value of an input, 
and measuring the color. It is equivalent to preparing a live-data pair. It is YMCK% of an input and the 
chromaticity coordinate of an output, for example, L* a* b*, by outputting a color patch in input 
[ YMCK% of ] various combination in the case of 4 color color printer, and measuring the color to a 
pan. It is equivalent to preparing a live-data pair. 

[0033] Next, at the forecast-ed preparation production process of SI 2, when the forecast-ed stated by the 
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above-mentioned (6) types and (7) formulas, i.e., a forecast-ed, is an input signal, (xlj, x2j, x3j) Or 
when a forecast-ed is an output signal about (xlj, x2j, x3j, x4j) again, only a required number prepares 
(ylj,y2j,y3j)or(ylj,y2j,y3j),and(x4j). . 
[0034] Next, the initial value of the coefficient Wij of weightmg of an n-tuple shown by (3) formulas is 
computed at the initial value calculation production process of the coefficient of weighting of SI 3. It is 
important for initial value to make it the value possible nearest to the value finally determined by the 
successive approximation method in order to bring convergence forward, for example, the input space 
which standardized the difference for every signal component of the live data of the forecast-ed input 
signal shown in (8) types and (9) types, and an input signal when a forecast-ed was an input signal - 
difference - weight Wlij of the monotonically decreasing function Fij of the Euclidean distance in the 
standardization input signal space for which it asks from a component It uses. The monotonically 
decreasing function Gij shown in (10) types and (1 1) types since the component of a matrix was not 
determined in this phase does not use, but is weight W 1 ij . It considers as initial value. 
[0035] Output space is substituted when a forecast-ed is an output signal, the output space which 
standardized the difference for every signal component with the live data of the forecast-ed output signal 
shown in (13) types shown below when an input was a three dimension -- difference ~ weight W2ij' of 
the monotonically decreasing function Gij of the Euclidean distance in the standardization input signal 
space for which it asks from a component is used, moreover -- an input -- four - a dimension - it is ~ a 
case - (-- 13 --) - a formula - having been shown -- weight - W - two - ij ~ ' - (-- 14 --) - a formula 
- having been shown - an input space - one - a ** ~ difference -- a component - having 
standardized - an input signal ~ space - it can set - distance -- a monotonically decreasing function - 
Gij - depending - weight - W - two - ij - 1 - (-- 
[0036] 

[Equation 5] „ 

W'^ = Gij(((y u - yij)/yu>? + - 2/ 2j )/2/2o) 

+((yzi - y 3j )/y3o) 2 ) (13) 



W' Ui = Fo(((x 4i -^i)Mo) 2 ) ( 14 ) 

W'm = H ( w w w *>) (15) 

[0037] Next, the square sum Ej of the Euclidean distance by which weighting of two or more output live 
data corresponding to it was carried out to the output forecast calculated at the component calculation 
production process of the matrix of SI 4 using the matrix from the live data of two or more input signals 
shown by (3) formulas The least square method is used for the basis of the conditions of making it mm, 
and each components ml 1 and ml2 of a matrix and ... are computed for a while. The coefficient Wij of 
weighting used here is the weight decided at the initial value calculation production process of the 
coefficient of weighting of S 13. . 
[0038] Next, a forecast is computed for a while from a forecast-ed at the forecast calculation production 
process of S15 using (6) types or (7) types. A forecast-ed An input signal (xlj, x2j, x3j) In (xlj, x2j, x3j, 
x4j), or (ylj, y2j, y3j) a forecast-ed - an output signal (ylj, y2j, y3j) -- with or (ylj, y2j, y3j), when it is 
(x4j), it asks for (xlj, x2j, x3j) using each components ml 1 and ml2 of the matrix searched for at the 
component calculation production process of the matrix of SI 4, and ... 

[0039] At the re-calculation production process of the coefficient of weighting of SI 6, next, (8) types or 
(9) types, the input space which standardized the difference for every signal component with the live 
data of the input signal of a forecast-ed or a forecast, and an input signal - difference - weight Wlij of 
the monotonically decreasing function Fij of the Euclidean distance in the standardization input signal 
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space for which it asks from a component (10) A formula or (1 1) types, weight W2ij of the 
monotonically decreasing function Gij of the Euclidean distance in the standardization output signal 
space for which it standardizes and asks after changing into the component of the output space which 
considered sensitivity using the component of the matrix by the difference for every signal component 
with the live data of the input signal of a forecast-ed or a forecast, and an input signal (12) -- the 
composite function H of a formula -- compounding -- the base of an n-tuple -- weight W12ij is again 
calculated to a data pair. 

[0040] Next, the square sum Ej of the Euclidean distance by which weighting of two or more output live 
data corresponding to it was carried out to the output forecast calculated at the component calculation 
production process of the matrix of SI 7 using the matrix from the live data of two or more input signals 
shown by (3) formulas The least square method is used for the basis of the conditions of making it min, 
and each components m 1 1 and m 1 2 of a matrix and ... are computed again. 

[0041] Next, a forecast is again computed from a forecast-ed at the forecast re-calculation production 
process of S18 using (6) types or (7) types, a forecast-ed ~ an input signal (xlj, x2j, x3j) - or (xlj, x2j, 
x3j, x4j) ~ it is ~ a case (ylj, y2j, y3j) a forecast-ed ~ an output signal (ylj, y2j, y3j) - with or (ylj, y2j, 
y3j), when it is (x4j), it asks for (xlj, x2j, x3j) using each components mil and ml2 of the matrix 
searched for at the component re-calculation production process of the matrix of S17, and ... 
[0042] Next, it judges whether at the convergence test production process of the forecast of SI 9, the 
forecast calculated at the forecast re-calculation production process of S18 was compared with the 
forecast calculated at the forecast calculation production process of S15, or the forecast re-calculation 
production process of SI 8 of 1 time ago, and it was completed by the forecast. As the method of a 
judgment, it asks for the difference between this forecast and a previous forecast for every component, 
and when the absolute value is smaller than the threshold for every component set up beforehand, it can 
judge with converging, for example. Or when the Euclidean distance between this forecast and a 
previous forecast is smaller than the threshold set up beforehand, it can also judge with converging. 
Here, if it judges with not progressing and converging on the prediction termination judging production 
process of S22 by making into a final forecast the forecast calculated at the forecast re-calculation 
production process of SI 8 this time if it judges with converging, it will progress to the re-calculation 
judging production process of the forecast of S20. 

[0043] The re-calculation judging production process of the forecast of S20 is a production process 
which judges whether the convergence direction of a forecast and the count of convergence count are 
investigated, and convergence count is made to perform further. When not emitting and converging 
according to this production process around extremal value with a forecast, he is trying not to lapse into 
the count loop of an infinity time. For example, when this production process is the 1st time, it judges 
with re-calculating a forecast unconditionally, and judges with 2nd henceforth re-calculating a forecast, 
when investigating the relation between a forecast and the forecast of 1 time ago, and the relation 
between the forecast of 1 time ago, and the forecast of 2 times ago for every component and having 
shifted in the same direction. What is necessary is just to judge with not re-calculating, since it does not 
emit and converge, when the count was counted when having shifted in the different direction, it judges 
with re-calculating when it is below the threshold that counted value set up beforehand, and a threshold 
is exceeded. Here, if it judges with re-calculating a forecast, it will progress to the re-calculation 
production process of the coefficient of weighting of SI 6, and if it judges that convergence count does 
not perform a repeat and a re-calculation, it will progress to the forecast terminal decision production 
process of S21. 

[0044] The forecast terminal decision production process of S21 is a production process which 
determines the value most appropriate for it as a final forecast, when not emitting and converging 
around extremal value with a forecast. The value most appropriate for it is determined as a final forecast 
- memorize the forecast which re-calculated, for example, take the average of the total forecast after 
beginning emission, or with a forecast and a forecast [ 1st ] difference chooses the smallest forecast. 
[0045] The following prediction termination judging production process of S22 is a production process a 
forecast judges whether calculation was completed or not about all forecasts-ed to be, in not ending, it 
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returns to the initial value calculation production process of the coefficient of weighting of SI 3, and, in a 
repeat and termination, all color transfer-characteristics prediction processings of a color picture I/O 
device end prediction about the following forecast-ed. 

[0046] Next, the gestalt of operation of the 2nd of the color transfer-characteristics prediction method of 
this invention is explained. With the gestalt of this 2nd operation, as a coefficient Wij of weighting of 
(3) types the input space which standardized the difference for every signal component with the live data 
of the input signal of the forecast-ed or forecast used with the gestalt of the 1st operation of a ****, and 
an input signal — difference - with the monotonically decreasing function Fij of the Euclidean distance 
in the standardization input signal space for which it asked from the component Besides the 
monotonically decreasing function Gij of the Euclidean distance in the standardization output signal 
space which was standardized after changing into the component of the output space which considered 
sensitivity using the component of this matrix by the difference for every signal component with the live 
data of the input signal of a forecast-ed or a forecast, and an input signal and for which it asked Data 
precision function Ji The example to be used is shown. This data precision function Ji They are the two 
above-mentioned kinds of weight Wlij(s) to the live data of an input signal. W2ij It is the function 
calculated from the interrelation of the output forecast used and predicted and the live data of an output 
signal, and the probability containing a noise of each live-data pair is expressed. 
[0047] the input space which standardized the difference for every signal component with the live data 
of the input signal of a forecast-ed or a forecast here - difference - with the monotonically decreasing 
function Fij of the Euclidean distance in the standardization input signal space for which it asked from 
the component The monotonically decreasing function Gij of the Euclidean distance in the 
standardization output signal space which was standardized after changing into the component of the 
output space which considered sensitivity using the component of this matrix by the difference for every 
signal component with the live data of the input signal of a forecast-ed or a forecast and an input signal 
and for which it asked data precision function Ji calculated from the interrelation of the output forecast 
predicted to the live data of an input signal using two kinds of such weight, and the live data of an 
output signal since it is the same as what was shown with the gestalt of the 1st operation of a **** The 
case where an output signal is a three dimension is taken for an example. ******** - It explains using 
(23) types from (16) types shown below. 
[0048] 
[Equation 6] 
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Aya = y' n - yn (16) 

&Vi2 = Vi2 ~ Vi2 ( 17 ) 

Ayo = 24-2*3 (18) 

= [EW«i-w)))/[E(H)) (19) 

= (EW(^-2/2 fc )))/(E(V fc )) (20) 

A»« = ^(^(^-^/^(V*)) (21) 



W si = Ji(In {Ay iu Ay ial ),I i2 {Ayi2, Aj/ ts2 ), Ix(&ya, Aj/^)) 

(22) 

^i23,i = H(H^,W 2li ,W 3l ) (23) 

[0049] (16) the differential factor of the output forecast in a live-data pair, and ****** data (18) types 
show the differential factor of the output forecast (yiT, and yi2' and yi3 - ') in the i-th live-data pair, and 
output live data (yil, yi2, yi3) from the formula, and excluding [ (19) types to (21) types ] the i-th live- 
data pair — weighting factor Vk What carried out the weighting average is shown, the input space where 
weight here standardized the difference for every signal component of the input signal of the i-th live 
data, and the input signal of live data other than the i-th ~ difference - the monotonically decreasing 
function of the Euclidean distance in the standardization input signal space for which it asks from a 
component — then, it is good. 

[0050] The following (22) types show the data precision function which the differential factor to (16) to 
(21) type is made to correspond for every differential factor, and asks for it. data precision function Ji 
shown here the difference of the output forecast and output live data about the i-th live-data pair — the 
difference of a weighting average with a direction, the output forecast about the live-data pair of a live- 
data pair except the live-data pair of magnitude **yil of the difference, **yi2, and 3 or i-th **yi, and 
output live data — magnitude **yisl of a direction and its difference - **yis2 and **yis3 from — it is 
the feature to use the functions Iil, Ii2, and Ii3 which show the probability of a live-data pair which took 
such a small value that the magnitude of a difference will become large about such a big value that the 
magnitude of a difference becomes large if both are the same directions if it is a different direction. 
[0051] The following (23) types are weight Wlij which used two monotonically decreasing functions 
Fij and Gij, and W2ij. Data precision function Ji Used weight W3i is compounded with a composite 
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function H, and it is coefficient W123ij of weighting. Creating is shown. Here, the synthetic method that 
the feature that two monotonically decreasing functions which are the element are monotonically 
decreasing functions of each Euclidean distance is not broken down, and the synthetic method that the 
feature that a data precision function shows the probability of a live-data pair is not broken down need to 
be used for the method of compounding, for example, a synthetic method like the sum or a product can 
be used for it. 

[0052] The flow chart which shows an example of actuation [ in / in drawing 2 / the gestalt of operation 
of the 2nd of the color transfer-characteristics prediction method of this invention ], and drawing 3 are 
flow charts which show the details of a data precision function preparation production process. Here the 
example which enabled it to use a data precision function alternatively is shown, and when not using a 
data precision function is chosen, it is the same as the gestalt of the 1st operation of a ****. Since the 
I/O live-data pair preparation production process of SI 1 is completely the same as that of the case where 
weight by the data precision function is not used, explanation is omitted. Next, it judges whether at the 
data precision function weight use judging production process of S23, the coefficient of weighting using 
a data precision function is used. In not using, it progresses to SI 2, and processing explained with the 
gestalt of the 1st operation of a **** is performed. In using the coefficient of weighting using a data 
precision function, it progresses to the data precision function weight preparation production process of 
S24. 

[0053] The data precision function weight preparation production process of S24 is shown in drawing 3 . 
Hereafter, it explains using drawing 3 . At the forecast-ed setting production process of S3 1, all the input 
live data prepared at the I/O live-data pair preparation production process of SI 1 of drawing 2 are set up 
as a forecast-ed, and the preparations which predict the output value over input live data less than [ S32 ] 
are made. 

[0054] Next, the initial value of the coefficient Wij of weighting of the n-tuple stated by (3) formulas is 
computed at the initial value decision production process of the coefficient of weighting of S32. the 
input space which standardized the difference for every signal component with the live data of the 
forecast-ed input signal which showed initial value in (8) types, and an input signal ~ difference — 
weight Wlij of the monotonically decreasing function Fij of the Euclidean distance in the 
standardization input signal space for which it asks from a component What is necessary is just to use. 
[0055] Next, the square sum Ej of the Euclidean distance by which weighting of two or more output live 
data corresponding to it was carried out to the output forecast calculated at the component calculation 
production process of the matrix of S3 3 using the matrix from the live data of two or more input signals 
shown by (3) formulas The least square method is used for the basis of the conditions of making it min, 
and each components ml 1 and ml2 of a matrix and ... are computed for a while. The coefficient Wij of 
weighting here is the coefficient Wij of weighting decided at the initial value calculation production 
process of the coefficient of weighting of S32. 

[0056] Next, a forecast is computed for a while from a forecast-ed at the forecast calculation production 
process of S34 using (6) types. Since a forecast-ed is an input signal (xlj, x2j, x3j), it asks for (ylj, y2j, 
y3j) using each components ml 1 and ml 2 of the matrix searched for at the component calculation 
production process of the matrix of S33, and ... 

[0057] At next, the re-calculation production process of the coefficient of weighting of S35 (8) - the 
input space which standardized the difference for every signal component with the live data of the input 
signal of the forecast-ed shown in a formula, or a forecast, and an input signal — difference — weight 
Wlij of the monotonically decreasing function Fij of the Euclidean distance in the standardization input 
signal space for which it asks from a component (10) It standardizes, after changing the difference for 
every signal component with the live data of the input signal of the forecast-ed shown in a formula, or a 
forecast, and an input signal into the component of the output space which considered sensitivity using 
the component of this matrix, weight W2ij of the monotonically decreasing function Gij of the 
Euclidean distance in the standardization output signal space for which it asks weight W12ij 
compounded with the composite function H shown in (12) types - the base of an n-tuple - it calculates 
again to a data pair. 
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[0058] Next, the square sum Ej of the Euclidean distance by which weighting of two or more output live 
data corresponding to it was carried out to the output forecast calculated at the component re-calculation 
production process of the matrix of S36 using the matrix from the live data of two or more input signals 
shown by (3) formulas The least square method is used for the basis of the conditions of making it min, 
and each components mil and ml 2 of a matrix and ... are computed again. 

[0059] Next, a forecast is again computed from a forecast-ed at the forecast re-calculation production 
process of S37 using (6) types. Since a forecast-ed is an input signal (xlj, x2j, x3j), it asks for (ylj, y2j, 
y3j) using each components ml 1 and ml2 of the matrix searched for at the component calculation 
production process of the matrix of S36, and ... 

[0060] Next, it judges whether at the convergence test production process of the forecast of S3 8, the 
forecast calculated at the forecast re-degree calculation production process of S3 7 was compared with 
the forecast calculated at the forecast calculation production process of S34, or the forecast re- 
calculation production process of S37 of 1 time ago, and it was completed by the forecast. As the 
method of a judgment, it asks for the difference between this forecast and the last forecast for every 
component, and when the absolute value is smaller than the threshold for every component set up 
beforehand, it can judge with converging, for example. Or when the Euclidean distance between this 
forecast and the last forecast is smaller than the threshold set up beforehand, it can also judge with 
converging. Here, if it judges with not progressing and converging on the prediction termination judging 
production process of S41 by making into a final forecast the forecast calculated at the forecast re- 
calculation production process of S37 if it judges with converging, it will progress to the re-calculation 
judging production process of the forecast of S39. 

[0061] The re-calculation judging production process of the forecast of S3 9 is a production process 
which judges whether the convergence direction of a forecast and the count of convergence count are 
investigated, and convergence count is made to perform further. When not emitting and converging by 
this processing around extremal value with a forecast, he is trying not to lapse into the count loop of an 
infinity time. For example, when this production process is the 1st time, it judges with re-calculating a 
forecast unconditionally, and judges with 2nd henceforth re-calculating a forecast, when investigating 
the relation between this forecast and the last forecast, and the relation between the last forecast and the 
forecast before last for every component and having shifted in the same direction. What is necessary is 
just to judge with not re-calculating, since it does not emit and converge, when the count was counted 
when having shifted in the different direction, it judges with re-calculating when it is below the 
threshold that counted value set up beforehand, and a threshold is exceeded. Here, if it judges with re- 
calculating a forecast, it will progress to the re-calculation production process of the weighting factor of 
S35, and if it judges that convergence count does not perform a repeat and a re-calculation, it will 
progress to the forecast terminal decision production process of S40. 

[0062] The forecast terminal decision production process of S40 is a production process which 
determines the value most appropriate for it as a final forecast, when not emitting and converging 
around extremal value with a forecast. The average of the total forecast after memorizing the forecast 
which re-calculated, for example, beginning is taken, or difference with this forecast determines the 
value most appropriate for it by choosing the smallest forecast etc. as a final forecast. 
[0063] It is the production process which judges whether calculation of a forecast ended the following 
prediction termination judging production process of S41 about all forecasts-ed, and in not ending, it 
returns to the initial value calculation production process of the weighting factor of S32, and, in a repeat 
and termination, prediction is progressed to the data precision function weight calculation production 
process of S42 about the following forecast-ed. 

[0064] The following data precision function weight calculation production process of S42 is the data 
precision function Ji shown in (22) types from (16) types. Data precision function weight W3i is 
computed by following, the difference of an output forecast and output live data concerning [ a data 
precision function ] an object live-data pair - the magnitude of a direction and its difference - the 
difference of a weighting average with the output forecast and output live data about the live-data pair of 
a live-data pair except an object live-data pair - from the magnitude of a direction and its difference It is 
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the function which took such a small value that the magnitude of a difference will become large about 
such a big value that the magnitude of a difference becomes large if both are the same directions if it is a 
different direction and which shows the probability of a live-data pair. 

[0065] It is the production process which judges whether the weight calculation by the data precision 
function ended the following data precision function weight calculation termination judging production 
process of S43 about all I/O live-data pairs, and in not ending, it returns to the data precision function 
weight calculation production process of S42, and it repeats the weight calculation by the data precision 
function about the following I/O live-data pair. In termination, it progresses to the forecast-ed 
preparation production process of S12 of drawing 2 . 

[0066] Since the prediction production processes from S12 to S22 of drawin g 2 are also almost the same 
as when [ also when using data precision function weight ] not using, only a different portion adds 
explanation. At the initial value calculation production process of the coefficient of weighting of SI 3, 
the initial value of the weighting factor Wij of an n-tuple shown by (3) formulas is computed, the input 
space where the initial value in this case standardized the difference for every signal component with the 
live data of the forecast-ed input signal shown in (8) types or (9) types, and an input signal when a 
forecast-ed was an input signal — difference ~ weight Wlij of the monotonically decreasing function Fij 
of the Euclidean distance in the standardization input signal space for which it asks from a component 
What compounded data precision function weight W3i by (24) formulas shown below is used, the 
output space which standardized the difference for every signal component with the live data of the 
forecast-ed output signal shown in (13) types, and an output signal when a forecast-ed was an output 
signal and an input was a three dimension - difference - what compounded weight W2ij' of the 
monotonically decreasing function Gij of the Euclidean distance in the standardization input signal 
space for which it asks from a component, and data precision function weight W3i by (25) formulas 
shown below is used, an input - four -- a dimension -- it is -- a case - (— 13 — ) - a formula - weight -- 
W — two - ij — ' - (-- 14 --) — a formula — having been shown - an input ~ space — one ~ a ** — 
difference - a component -- having standardized ~ an input signal ~ space -- it can set - distance - a 
monotonically decreasing function - Fij - depending -- weight ~ W - one ~ ij - ' - data -- precision - 
a function — weight 
[0067] 
[Equation 7] 

W[ 3ij = H(W Uj ,W 5i ) (24) 
WUij = H(W^W 3i ) (25) 
W[ 23ij = H(W[ ijy W^W 3i ) (26) 

[0068] At the re-calculation production process of the coefficient of weighting of S16 (8) - the input 
space which standardized the difference for every signal component with the live data of the input signal 
of the forecast-ed shown in a formula or (9) types, or a forecast, and an input signal - difference - 
weight Wlij by the monotonically decreasing function Fij of the Euclidean distance in the 
standardization input signal space for which it asks from a component (10) It standardizes, after 
changing the difference for every signal component with the live data of the input signal of the forecast- 
ed shown in a formula or (1 1) types, or a forecast, and an input signal into the component of the output 
space which considered sensitivity using the component of this matrix. Weight W2ij by the 
monotonically decreasing function Gij of the Euclidean distance in the standardization output signal 
space for which it asks Data precision function Ji shown in (22) types Weight W3i to depend is 
compounded with the composite function H of (23) types, and it is weight W123ij to the live-data pair of 
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an n-tuple. It calculates again. 

[0069] It is also completely the same as when [ also when using data precision function weight ] not 
using, and prediction in the case of using data precision function weight can be carried out by passing 
through the prediction production processes from S12 to S22 except SI 3 and S 16 of a more than. 
[0070] Next, the gestalt of operation of the 3rd of the color transfer-characteristics prediction method of 
this invention is explained, the input space standardized by breaking the difference for every signal 
component with the live data of the input signal of a forecast-ed, and an input signal used with the 
gestalt of the 1st operation as a coefficient Wij of weighting of (3) types by the constant for every signal 
component by the gestalt of this 3rd operation — difference — the example only using the monotonically 
decreasing function Fij of the Euclidean distance in the standardization input signal space for which it 
asks from a component is shown. 

[0071] here - -ed -- a forecast -- an input signal - an input signal — live data - a signal component -- 
every -- a difference -- each — a component — every - that - a constant — dividing - things — having 
standardized - an input -- space - difference -- a component -- from - asking -- standardization - an 
input signal - space ~ it can set — Euclidean distance ~ a monotonically decreasing function — Fij -- a 
**** - the -- one - operation - a gestalt — having been shown -- (— eight — ) — a formula — or Since a 
monotonically decreasing function Fij is dependent on the input signal which is a forecast when 
searching for the input signal which is a forecast corresponding to the output signal of the arbitration 
which is a forecast-ed, the color transfer-characteristics prediction method becomes completely the same 
as the method shown in drawing 1 as a coefficient Wij of weighting of (3) types except using Wij shown 
in (8) types or (9) types instead of using W12ij of (12) types. Therefore, the case where the output signal 
which is a forecast corresponding to the input signal of the arbitration which is a forecast-ed here is 
searched for is explained. 

[0072] Drawing 4 is a flow chart which shows an example of the actuation in the case of searching for 
the output signal which is a forecast corresponding to the input signal of the arbitration which is a 
forecast-ed in the gestalt of operation of the 3rd of the color transfer-characteristics prediction method of 
this invention. Although it is the I/O live-data pair preparation production process of S51, since this is 
completely the same as that of the I/O live-data pair preparation production process of SI 1 shown in 
drawin g 1 at first, explanation is omitted. 

[0073] Next, at the forecast-ed preparation production process of S52, since it is the case where a 
forecast-ed is an input signal, only a required number prepares the input signal (xlj, x2j, x3j) or input 
signal (xlj, x2j, x3j, x4j) of a forecast-ed stated by above-mentioned (6) types and (7) formulas. 
[0074] Next, the coefficient Wij of weighting of the n-tuple stated by (3) formulas is computed at the 
calculation production process of the weighting coefficient of S 5 3. the input space which is the case 
where a forecast-ed is an input signal, and was standardized as a coefficient of weighting by dividing the 
difference for every signal component with the live data of the input signal of the forecast-ed shown in 
(8) types or (9) types, and an input signal by the constant for every signal component — difference - it 
can ask uniquely that what is necessary is to ask only for the monotonically decreasing function Fij of 
the Euclidean distance in the standardization input signal space for which it asks from a component. 
[0075] Next, the square sum Ej of the Euclidean distance to which weighting of the live data of two or 
more output signals corresponding to the output forecast and it for which it asked at the component 
calculation production process of the matrix of S54 using the matrix from the live data of two or more 
input signals shown by (3) formulas was carried out On the basis of the conditions of making it min, 
each components mil and ml 2 of a matrix and ... are computed using the well-known method of the 
least square method. 

[0076] Next, a forecast is computed from a forecast-ed at the forecast calculation production process of 
S55 using (6) types or (7) types. Since it is the case where a forecast-ed is an input signal (xlj, x2j, x3j) 
or an input signal (xlj, x2j, x3j, x4j), an output signal (ylj, y2j, y3j) is searched for using each 
components ml 1 and ml 2 of the matrix searched for at the component calculation production process of 
the matrix of S54, and ... 

[0077] The following prediction termination judging production process of S56 is a production process a 
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forecast judges whether calculation was completed or not about all forecasts-ed to be, in not ending, it 
returns to the calculation production process of the weighting coefficient of S53, and it repeats 
prediction about the following forecast-ed. All color transfer-characteristics prediction of the color 
picture I/O device in the case of using only a monotonically decreasing function Fij in termination is 
ended. 

[0078] Since the coefficient of weighting can be uniquely defined when using only a monotonically 

decreasing function Fij by the case where the output signal which is a forecast corresponding to the 

input signal of the arbitration which is a forecast-ed is searched for so that it may understand by the 

above explanation, predicting without using the successive approximation method is possible. 

[0079] Next, the gestalt of operation of the 4th of the color transfer-characteristics prediction method of 

this invention is explained. With the gestalt of this 4th operation, as a coefficient Wij of weighting of (3) 

types the input space standardized by dividing the difference for every signal component with the live 

data of the input signal of a forecast-ed, and an input signal used with the gestalt of the 1st operation by 

the constant for every signal component - difference - with the monotonically decreasing function Fij 

of the Euclidean distance in the standardization input signal space for which it asks from a component 

Data precision function Ji used with the gestalt of the 2nd operation The example to be used is shown. 

[0080] the input space standardized here by dividing the difference for every signal component with the 

live data of the input signal of a forecast-ed, and an input signal by the constant for every signal 

component — difference - the monotonically decreasing function Fij of the Euclidean distance in the 

standardization input signal space for which it asks from a component It is the same as (8) types or (9) 

types which were shown with the gestalt of the 1st operation of a ****, and is the data precision 

function Ji. It is the same as the (16) type - (22) type shown with the gestalt of the 2nd operation of a 
**** 

[0081] It is Wlij shown in (8) types or (9) types instead of W12ij of (12) types being used for the color 
transfer-characteristics prediction method as a coefficient Wij of weighting of (3) types since the 
monotonically decreasing function Fij was dependent on the input signal which is a forecast when the 
input signal which is a forecast corresponding to the output signal of the arbitration which is a forecast- 
ed was searched for. Except using, it becomes completely the same as the method shown in drawing 2 or 
drawing 3 . Therefore, the case where the output signal which is a forecast corresponding to the input 
signal of the arbitration which is a forecast-ed here is searched for is explained. 

[0082] The flow chart and drawin g 6 which show an example of actuation in case drawing 5 searches for 
the output signal which is a forecast corresponding to the input signal of the arbitration which is a 
forecast-ed in the gestalt of operation of the 4th of the color transfer-characteristics prediction method of 
this invention are a flow chart which similarly shows the details of a data precision function preparation 
production process. Here, the example which enabled it to use a data precision function alternatively is 
shown, and when not using a data precision function is chosen, it becomes completely the same as the 
gestalt of the 3rd operation of a ****. Although it is the I/O live-data pair preparation production 
process of S61 at first, since this is completely the same as that of the I/O live-data pair preparation 
production process of SI 1 in drawing 1 , explanation is omitted. 

[0083] Next, it judges whether at the data precision function weight use judging production process of 
S62, the weighting coefficient using a data precision function is used. When not using, it progresses to 
S64, and processing explained with the gestalt of the 3rd operation of a **** is performed. In using the 
coefficient of weighting using a data precision function, it progresses to the data precision function 
weight preparation production process of S63. Since the data precision function weight preparation 
production process of S63 is shown in drawing 6 , it explains using drawing 6 below. 
[0084] At the forecast-ed setting production process of S71, all the I/O live data prepared at the I/O live- 
data pair preparation production process of S61 of drawing 5 are set up as a forecast-ed, and the 
preparations for predicting the output value over I/O live data less than [ S72 ] are made. 
[0085] Next, at the decision production process of the coefficient of weighting of S72, the coefficient 
Wij of weighting of the n-tuple stated by (3) formulas is computed, the input space which is the case 
where a forecast-ed is an input signal, and was standardized as a coefficient of weighting by dividing the 
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difference for every signal component with the live data of the forecast-ed input signal shown in (8) 
types or (9) types, and an input signal by the constant for every signal component - difference — it can 
set uniquely that what is necessary is to ask only for the monotonically decreasing function Fij of the 
Euclidean distance in the standardization input signal space for which it asks from a component. 
[0086] Next, the square sum Ej of the Euclidean distance by which weighting of the live data of two or 
more output signals corresponding to it was carried out to the output forecast at the component 
calculation production process of the matrix of S73 using the matrix from the live data of two or more 
input signals shown by (3) formulas On the basis of the conditions of making it min, the components 
ml 1 and ml2 of a matrix and ... are computed using the well-known method of the least square method. 
The coefficient of weighting here is a coefficient of weighting decided at the calculation production 
process of the coefficient of weighting of S72. 

[0087] Next, at the forecast calculation production process of S74, a forecast is computed from a 
forecast-ed using (6) types or (7) types. Since a forecast-ed is an input signal (xlj, x2j, x3j), an output 
signal (ylj, y2j, y3j) is searched for using each components ml 1 and ml2 of the matrix searched for at 
the component calculation production process of the matrix of S73, and ... 

[0088] The following prediction termination judging production process of S75 is a production process a 
forecast judges whether calculation was completed or not about all forecasts-ed to be. In not ending, it 
returns to the calculation production process of the coefficient of weighting of S72, and it repeats 
prediction about the following forecast-ed. In termination, it progresses to the data precision function 
weight calculation production process of S76. 

[0089] The data precision function weight calculation production process of S76 computes data 
precision function weight according to the data precision function shown in (19) types from (16) types, 
the difference of the output forecast concerning [ a data precision function ] an object live-data pair and 
the live data of an output signal - the magnitude of a direction and its difference - the difference of a 
weighting average with the output forecast about the live-data pair of a live-data pair and the live data of 
an output signal except an object live-data pair ~ from the magnitude of a direction and its difference It 
is the function which was made to make such a big value that the magnitude of a difference become 
large if both are the same directions a value smaller if it is a different direction, as the magnitude of a 
difference will become large and which shows the probability of a live-data pair. 
[0090] The following data precision function weight calculation termination production process of S77 
is a production process which judges whether the weight calculation by the data precision function was 
completed about all I/O live-data pairs. In not ending, it returns to the data precision function weight 
calculation production process of S76, and, in a repeat and termination, calculation of data precision 
function weight is progressed to the forecast-ed preparation production process of S64 in drawing 5 
about the following I/O live-data pair. 

[0091] Since it is almost the same and only S64 differs also when data precision function weight is used 
for the prediction production processes from S64 to S68 of drawing 5 , and also when not using, only 
S64 adds explanation. 

[0092] At the calculation production process of the weighting factor of S64, the coefficient Wij of 
weighting of the n-tuple stated by (3) formulas is computed, in this case, the input space which 
standardized the difference for every signal component with the live data of the input signal which is the 
forecast-ed shown in (8) types or (9) types since a forecast-ed is an input signal, and an input signal - 
difference - it can set uniquely that what is necessary is just to use what compounded the weight of the 
monotonically decreasing function of Euclidean distance and data precision function weight in the 
standardization input signal space for which it asks from a component by (24) formulas. 
[0093] As mentioned above, except S64 of drawing 5 , it is also completely the same as when using data 
precision function weight and not using, and prediction using data precision function weight can be 
performed by passing through the prediction production processes from S61 to S68. Since the 
coefficient of weighting can be uniquely defined when using only a monotonically decreasing function 
Fij by the case where the output signal which is a forecast corresponding to the input signal of the 
arbitration which is a forecast-ed is searched for, even if it is the case where data precision function 
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weight is used, as it understands by the above explanation, predicting without using the successive 
approximation method is possible. 

[0094] As mentioned above, although four gestalten of operation of this invention were explained, in 
addition although the combination of the color transfer-characteristics prediction method shown in the 
data precision function shown in drawing 6 and drawing 2 is also possible, since it is the combination 
which can be guessed easily, explanation is omitted. 
[0095] 

[Example] Drawin g 7 is the block diagram showing the 1st example of the system which applies the 
color transfer-characteristics prediction method of this invention. As for a color scanner, and 102-104, 
101 are [ a 1 -dimensional table and 105 ] the matrix masking sections among drawing. Here, how to be 
on count and to check the precision of the method of determining the coefficient of color processing 
using the color transfer-characteristics prediction method of this invention, and its processing is 
explained, using a color scanner 101 as a color picture input unit. 

[0096] In the system shown in drawin g 7 , a color scanner 101 reads the given color picture, and outputs 
the signal of a RGB color space. As for the signal of the RGB color space outputted from a color 
scanner 101, in the 1 -dimensional tables 102-104, color adjustment is performed for every signal 
component. The purpose of the 1 -dimensional tables 102-104 is data of the RGB color space usually 
proportional to a reflection factor at the following matrix masking section 105 L* a* b* So that 
conversion precision may become good, when changing into the data of a color space In order to amend 
the difference in the sensitivity of each sensor of changing into data R'G'B' of another unit which is not a 
reflection factor, for example, units, such as concentration and lightness, and RGB of a color scanner 
101, a color scanner — 101 - gray — inputting - having had - the time — R — ' — G — 1 — B — ' — being 
the same ~ a value ~ becoming — as — conversion — giving -- things — it is . uniform-color-space L* 
which is not dependent on equipment in the matrix masking section 105 after changing 1 -dimensional 
table 102-R, G, and B each component of every by 104 - it changes and outputs to 'a^b*'. Uniform- 
color-space L* a* b* which measured the color picture made to read with a color scanner 101 with the 
colorimetry vessel etc. at this time Live data and the data of L* f a* f b* ? color space of the output of the 
matrix masking section 105 are made in agreement. The color picture data which can reproduce 
faithfully the color picture given to the color scanner 101 by this will be inputted. 
[0097] Drawing 8 is a flow chart which shows an example of the color processing coefficient decision in 
the 1st example of the system which applies the color transfer-characteristics prediction method of this 
invention, and check processing. Here, the prediction method which does not use a data precision 
function and which was explained with the gestalt of the 1st operation shall be used as a prediction 
process of color transfer characteristics. 

[0098] First, a color is beforehand measured at the I/O live-data pair preparation production process of 
S81, and it is uniform-color-space L* a* b*. A color scanner 101 is made to read the various color 
patches which obtained the data which can be set, and the signal in the RGB color space outputted from 
a color scanner 101 is measured. Thereby, it is uniform-color-space L*a* b* as live data of an input 
signal. The live-data pair of data and the live data in the RGB color space of the output of a color 
scanner 101 can be prepared. As for the color of various color patches, being uniformly distributed in 
the color space of an input is desirable, and the number usually ranks 100 to 1000th, although it is 
dependent on the precision of the prediction to desire. 

[0099] Since the noise according to causes, such as heterogeneity within the field of a color scanner and 
a time unstable property, in an output signal the time of reading of a color picture has ridden, it is good 
to devise so that these noises may be easy to be removed by the smoothing feature which the color 
property prediction method has. Moreover, the effect of a time unstable property is also good for a RGB 
value by making the array random at the time of reading of a color patch, making it the effect of the 
heterogeneity within a field ride on the data of a RGB color space at random, or dividing into a many 
times scan and extracting the data of a RGB color space to make it ride at random. Furthermore, after 
also taking the heterogeneity of the color patch itself into consideration, not being made to make small 
magnitude of the aperture at the time of measurement of a color too much or measuring the data of a 
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RGB color space, it is still better to devise in quest of the average of the field corresponding to the 
magnitude of aperture etc. 

[0100] Next, the configuration and parameter of the function of weighting at the time of color transfer- 
characteristics prediction are set up at the weighting function configuration / parameter setup production 
process of S82. For example, it sets up as a (27) - (29) type showing the function of weighting 
equivalent to (8), (10), and (12) types below. Here, for (Lj, aj, and bj), the live data of an input signal, 
and (L0, aO and bO) are [ the input signal of a forecast-ed or a forecast, and (Li, ai and bi) ] the constant 
of standardization of input signal space, the constant of standardization of (R0, GO, and B0) of output 
signal space, and the constant p decides a function configuration to be. 



W 2ij = lim^Li^L^f + im^a^a^f 

+(m l3 (b i - bj)) 2 )/(R,) 2 
+((m 21 (Li - Lj)) 2 + (m 22 (ai - aj )) 2 

+{m2z{bi ~ 6;)) 2 )/(Go) 2 
+((m 3 i(Li - Lj)) 2 + (m 32 0,' - a,)) 2 

+(m 33 (fr,- - bj)) 2 )/(B 0 ) 2 y + 1) (28) 

W l2ij = W Uj + W 2ij (29) 



[0102] The parameter of weighting means (L0, aO, bO), (R0, GO and B0), and three kinds of constants of 
p here, and the degree of a smoothing feature is adjusted. As it understands in the configuration of a 
formula, a smoothing feature is a parameter with which (L0, aO, bO), and (R0, GO and B0) become weak, 
so that, as for p, a value becomes large by a smoothing feature becoming strong, so that a value becomes 
large. Before these parameters perform prediction, they need to be determined by the degree of the noise 
of an I/O live-data pair. What is necessary is to strengthen a smoothing feature generally, if a noise is 
large, and just to weaken a smoothing feature, if a noise is small. However, since the data which usually 
judges the degree of a noise does not have only an I/O live-data pair and the true value is not necessarily 
known, a certain degree will be set up experientially. 

[0103] For example, the group of a certain parameter is assumed, the output to input live data is 
predicted, the difference of a forecast and output live data is seen, and there is the method of choosing 
the group of a certain suitable parameter. In this case, the difference of a forecast and output live data is 
that by which the error of the anticipation method itself and the error by the noise were compounded, 
and should just ask for the relation between this difference and the group of an optimal parameter by the 
simulation etc. beforehand. 

[0104] As another example, it is L* a* b* as a forecast-ed. It is L* in the gradation data of a color space, 
a * =b* =0 [ for example, ]. A forecast-ed which changes little by little is prepared, the group of a certain 
parameter is assumed, and the data of a RGB color space is predicted. And L* The data of the RGB 
color space of a prediction result to receive is graph-ized, and it judges whether the group of the 
assumed parameter is sufficient from the configuration. For example, a group of a parameter which the 
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result follows a live-data pair too much, and gives stronger smoothing since the smoothing feature is too 
weak when there is much unnatural extremal value is assumed again, and it judges similarly. The group 
of the optimal parameter can also be set up by repeating this several times. 

[0105] Next, two or more forecasts-ed for deciding the coefficient of color processing at the forecast-ed 
preparation production process of S83 are set up. Although it is dependent on the decision method of the 
color processing ARUKO rhythm to be used and its coefficient, a setup of this forecast-ed shall 
determine previously the coefficient of the 1 -dimensional tables 102-104 here, and then shall determine 
the coefficient of the matrix masking section 105. In this case, in order to decide the coefficient of the 1- 
dimensional tables 102-104, it is L* in a gray scale, a* =b* =0 [ i.e., ]. The forecast-ed which is 
changing equally from zero to 100 is set up. It will be satisfactory, if width of face of the step is made 
about into 100/256 when using a 8-bit table. Moreover, if it is the premise of wanting to make the color 
difference small equally in a color space in order to decide the coefficient of the matrix masking section 
105, it is L* a* b*. Lattice point data, for example, the lattice point data of ten gaps, is set up as a 
forecast-ed. The color may be added as long as there are colors of wanting to make especially the color 
difference small, such as people's complexion. 

[0106] Next, prediction is performed at the forecast calculation production process of S84 about two or 
more forecasts-ed set up at the forecast-ed preparation production process of S83. Since it is completely 
the same as that of the method explained by drawing 1 , the method is omitted. The gray scale and L* 
a*b* which are a forecast-ed by this The forecast in the RGB color space to lattice point data is 
obtained. 

[0107] Next, a data pair required in order to decide the coefficient of color processing is elected from the 
data pair of two or more forecasts-ed and a forecast obtained at the production process to S84 at the data 
election production process for color processing coefficient decision of S85. The criteria of election are 
whether the close forecast of a RGB color space is in the range treated by color processing of this 
system. That is, the data of a RGB color space is usually a gray scale and L* a* b*, although it can 
express with the unit of a reflection factor. The data which the forecast of the RGB color space of lattice 
point data has for less than 0 and 100% may be contained. What is necessary is to remove those data and 
just to elect the data pair to which the close forecast of a RGB color space is in 0 to 100% of range, since 
those data is unnecessary as data for color processing coefficient decision. 

[0108] Next, the coefficient for actual color processing is determined at the color processing coefficient 
decision production process of S86. First, the coefficient of the 1-dimensional tables 102-104 is 
determined using the data of a gray scale. Here, the 1-dimensional tables 102-104 shall change the data 
of the RGB color space proportional to a reflection factor into the data of the R f G'B f color space which 
makes lightness an unit. The coefficient of the 1-dimensional tables 102-104 corresponding to each 
signal component of R, G, and B is the value of the forecast of R, or G or B to an axis of ordinate in a 
horizontal axis L* A value is taken, the data pair of gray scale is plotted, the transform function from a 
RGB color space to a R'G'B 1 color space can be defined by approximating with the polygonal line 
between plots, and it can be determined by quantizing this transform function. 

[0109] Next, it is the coefficient of the matrix masking section 105 L* a* b* It determines using lattice 
point data. L* a* b* The forecast of the RGB color space corresponding to lattice point data It changes 
into the data of a R'G'B' color space using the 1-dimensional tables 102-104 decided previously. L* a* 
b* of the color picture inputted by considering the data of the R'G'B' color space as the input of the 
matrix masking section 105 By calculating the data of a color space recursively as an aim of an output 
using the least square method The coefficient of the matrix masking section 105 can be determined. 
[0110] At the last, the precision of color processing is checked at the color processing check production 
process of S87. L* a* b* of a color to check checking the precision of color processing It means 
investigating what has happened to the data of L*'a*'b*' color space which is an output after color 
processing to the data of a color space. Therefore, color L* a* b* of arbitration to check Consider as a 
forecast-ed and the data of the RGB color space which is the output of a color scanner 101 is predicted. 
It is changed into the data of a R'G'B' color space using the 1-dimensional tables 102-104, and it changes 
into the data of an L*'a*'b*' color space by the matrix masking section 105 further, and is L* a* b*. 
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What is necessary is just to search for the color difference with the data in a color space. It is possible to 
check the data of the RGB color space of arbitration for the same thing as a starting point. At this time, it 
is L* a* b* of the color picture given to a color scanner 101 by making the data of the RGB color space 
of arbitration into a forecast-ed. What is necessary is to predict the data of a color space and just to 
search for the color difference of what changed the data of a RGB color space into L**a*' b*' using the 
1 -dimensional tables 102-104 and the matrix king section 105. 

[0111] Drawin g 9 is the block diagram showing the 2nd example of the system which applies the color 
transfer-characteristics prediction method of this invention. As for a three-dimension table, and 112-115, 
1 1 1 are [ a 1 -dimensional table and 116] color printers among drawing. Here, how to be on count and to 
check the precision of the method of determining the coefficient of color processing using the color 
transfer-characteristics prediction method of this invention, and its processing is explained, using a color 
printer 1 16 as a color picture output unit. 

[0112] In the system shown in drawing 9 , the data of a color picture which should be recorded is given 
to the three-dimension table 1 1 1 as data of L*'a*'b*' color space. The three-dimension table 1 1 1 changes 
the data of an L*'a*'b*' color space into the data of the 4-dimensional color space of Y'M'C'K'. In that 
case, it asks by interpolation processing about the color which is not in a table. L* a* b* of the color 
picture recorded by the data of I^'a^b*' color space inputted, and the color printer 1 16 at this time The 
data in a color space changes so that it may be in agreement if possible. 

[0113] The data of the Y'M'C'K 1 color space outputted from the three-dimension table 1 1 1 is inputted 
into the 1 -dimensional tables 1 12-1 15 for every component. The 1 -dimensional each tables 1 12-1 15 
adjust each component, and output it to a color printer 1 16 as data of a YMCK color space. The purpose 
of these 1 -dimensional tables 1 12-1 15 is amending the local nonlinearity of the monochrome scale 
which cannot cope with it on the three-dimension table 111. 

[01 14] A color printer 116 records the data of the YMCK color space outputted from the 1 -dimensional 
each tables 1 12-1 15 on reception, and records a color picture on a record medium-ed. L* a* b* of the 
color picture recorded at this time The color picture which reproduced the given color picture data 
faithfully can be obtained by constituting the coefficient of the three-dimension table 111, and the 1- 
dimensional table 1 12 so that the value in a color space and the value in L*'a*'b*' color space of the 
color picture data inputted into the three-dimension table 111 may be in agreement. 
[0115] Drawing 10 is a flow chart which shows an example of the color processing coefficient decision 
in the 2nd example of the system which applies the color transfer-characteristics prediction method of 
this invention, and check processing. Here, as the gestalt of the 2nd operation of a **** explained as a 
prediction process, the prediction method which used the data precision function shall be used. 
[0116] First, it is L* a* b* of the data of the YMCK color space of an input side, and an output side by 
outputting a color patch in the various combination of the data of a YMCK color space, and measuring 
the color at the I/O live-data pair preparation production process of S91. A live-data pair with the data of 
a color space is prepared. Like the case of an above-mentioned color picture input unit, as for the color 
of various color patches, it is desirable to be uniformly distributed in the color space of an input side, for 
example, it is ideal to output a color patch in all the combination in every 10% of YMCK etc. However, 
since a number becomes large too much, even if it thins out using a rectangular table or the data of a 
YMCK color space changes, it is L* a* b* of an output side. By thinning out a field where the data of a 
color space seldom changes, it is practical to output about several 100 to several 1000 color patch. On 
the contrary, a portion from which the color transfer characteristics of a color printer 1 16 are changing 
rapidly may output a color patch at a finer gap. Moreover, it is good to also prepare the device which 
devises or copes with the heterogeneity of the color patch itself so that these noises may be easy to be 
removed by the smoothing feature which the color property prediction method has in the case of the 
output of a color patch, since the noise by the cause of the heterogeneity within the field of a printer and 
time instability has ridden on the output signal. Since it is the same as that of the case of an above- 
mentioned color picture input unit, those methods are omitted. 

[0117] Next, at the weighting function configuration / parameter setup production process of S92, the 
configuration and parameter of the function of weighting at the time of color transfer-characteristics 
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prediction are set up. For example, it sets up like (30) - (39) type showing the weighting function 
equivalent to (9), (1 1), and (16) - (23) type below. (Yj, Mj, Cj, and Kj) here The input signal of a 
forecast-ed or a forecast, The constant of standardization [ (Yi, Mi, Ci, and Ki) ] of the live data of an 
input signal, and (Y0, MO, CO and K0) of input signal space, The output signal of a forecast-ed or a 
forecast, and (Li, ai and bi) are the constants (Lj, aj, and bj) decide the live data of an output signal, and 
the constant of standardization of (L0, aO, and bO) of output signal space, and p decides function 
configurations to be. 
[0118] 
[Equation 9] 

W Uj = 1/(((W - *j)/*o) 2 + ((M - Mj)/M o y 

- C^/Co) 2 + ((*,- - K^flUff + 1) (30) 

w 2ij = i/((((m 11 (y;-y J -)) 2 + (mi 2 (M i -M i )) 2 

+(m 13 (C i - C^f + (muiKi - A',)) 2 )/^) 2 
+((m 21 (^ - lj)) 2 + (m 22 (M, - M } )f 
+(m 23 (C i - Cj)) 2 + (m 24 (/^ - ^)) 2 )/(a 0 ) 2 
+((m 3 i(li - Yj)) 2 + (m 32 (Mi - Mj)) 2 
+(m 33 (C,- - Cj)f + (m^Ki - Kj)) 2 )/(b Q ) 2 y + 1) 

(31) 

[Equation 10] 
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AL, = L'i-Li 
A a,- = a- - ai 
= b't-bi 



(32) 
(33) 
(34) 



AL is = 




(35) 



Aa ia = 




(36) 




(37) 



W 3i = ((\(AL ig + ALi) 2 x AL IS x AL,| l/4 /L 0 
+|(Aa,- a + Aa,-) 2 x Aa„ x Aa, | 1/4 /a 0 
+|(A6 ia + A& £ ) 2 x Ab ia x Ab^/boY + 1) 
/((\(AL ia - AX t ) 2 x AL £a x AL.f/VLo 
+|(Aa ifi - Aa,) 2 x Aa is x Aa,| 1/4 /a 0 
+|(A6,- 3 - A6,-) 2 x Ab is x A6 i | 1 / 4 /6 0 ) p + 1) (38) 

Wl23,; = ^X^X^i (39) 



[01 19] The parameter of weighting means (Y0, MO, CO, KO), (LO, aO and bO), and three kinds of 
constants of p, and the degree of a smoothing feature is adjusted here. As it understands in the 
configuration of a formula, a smoothing feature is a parameter with which (Y0, MO, CO, KO), and (LO, 
aO and bO) become weak, so that, as for p, a value becomes large by a smoothing feature becoming 
strong, so that a value becomes large. Before performing prediction, it is determined by the degree of the 
noise of an I/O live-data pair, but since it is the same as that of the case of an above-mentioned color 
picture input unit, these parameters omit the method. 

[0120] Next, the weight which is equivalent to (16) - (23) type at the data precision function weight 
preparation production process of S93, i.e., (32) - (39) type, to weight W123ij It calculates. Since it is 
completely the same as that of the method already explained in drawing 3 , the method is omitted here. 
The data precision function weight to the live-data pair of all I/O is determined by this production 
process. 
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[0121] Next, two or more forecasts-ed for deciding the coefficient of color processing at the forecast-ed 
preparation production process of S94 are set up. Although it is dependent on the decision method of the 
color processing ARUKO rhythm to be used and its coefficient, a setup of this forecast-ed shall 
determine previously the coefficient of the 1 -dimensional tables 112-115, next shall determine the 
coefficient of the three-dimension look-up table 1 1 1 here. In this case, in order to decide the coefficient 
of the 1 -dimensional tables 112-115 A monochrome scale, i.e., the thing from which Y is changing 
equally from 0% to 100% at M=C=K=0%, That from which K is changing equally from 0% to 100% at 
that from which C is changing equally from 0% to 100% at that from which M is changing equally from 
0% to 100% at Y=C=K=0%, and Y=M=K=0%, and Y=M=C=0% is set up as a forecast-ed. It will be 
satisfactory, if width of face of the step is made about into 100/256 when using a 8-bit table. 
[0122] Moreover, L* a* b* corresponding to the lattice point of the table in order to decide the 
coefficient of the three-dimension table 111 % value of the data of a color space and K corresponding to 
the color is set up as a forecast-ed. Various methods, such as a method of using the maximum ****, are 
devised, and the setting method of % value of this K is omitted here, although it is also possible to 
determine this using the prediction method of this invention. 

[0123] Next, prediction is performed at the forecast calculation production process of S95 about two or 
more forecasts-ed set up at the forecast-ed preparation production process of S94. The method is L* a* 
b* to the monochrome scale which is a forecast-ed although it omits since it is completely the same as 
that of the method explained by drawing 2 . The forecast of a color space, and L* a* b* of the lattice 
point The forecast of the YMC color space to the data of a color space and % data of K is obtained. 
[0124] Next, the coefficient of an actual color processor is determined at the color processing coefficient 
decision production process of S96. First, the coefficient of the 1 -dimensional tables 1 12-1 15 is 
determined using the data of a monochrome scale. The coefficient of the 1 -dimensional tables 112-115 
corresponding to each component of Y, M, C, and K L* a* b* which is the forecast of each 
monochrome scale at a horizontal axis The color difference value over 0% of chromaticity computed 
from the value of space % value of Y, M, or C or K is taken along an axis of ordinate, the data pair of a 
monochrome scale is plotted, by approximating with the polygonal line, the transform function from a 
Y'M'CK' color space to a YMCK color space can be defined between plots, and it can be determined by 
quantizing this transform function. 

[0125] Next, L* a* b* corresponding to the lattice point of the three-dimension table 111 for the 
coefficient of the three-dimension table 1 1 1 It determines using the data of a color space. For that 
purpose, what is necessary is to change into the data of a YTVTC'K' color space at reverse the 1- 
dimensional tables 112-115 which determined previously % value of K already determined as % value 
of the data corresponding to the lattice point of the three-dimension table 111, and the YMC color space 
for which it asked by prediction by solution Lycium chinense, and just to let the value be the coefficient 
of the corresponding lattice point. 

[0126] However, this method is the fundamental method of using, when the data of L*'a*'b*' color space 
of the lattice point of the three-dimension table 111 is in the color gamut which can reproduce a color 
printer 116. When the data besides that color gamut may be inputted into this system, it is necessary to 
extend this method. For example, the 1 -dimensional tables 1 12-1 15 are not created among 0 - 100%. If 
it creates temporarily in the field which covers all the values of L* f a* f b*' of the lattice point of the 
three-dimension table 1 1 1 using extrapolation and the same way of determining is adopted using this 
The data of a Y'M'CK'Y'M'CK including outside of that to which some data of color space exceeds less 
than 0% or 100%, i.e., color gamut, 1 color space can be found. Although this value can be made into the 
coefficient of the final lattice point after that combining the method of storing to 0 - 100%, the technique 
usually called Gamut compression, it omits for details. 

[0127] The precision of color processing is checked at the color processing coefficient check production 
process of S97 at the last. L* a* b* outputted from a color printer 1 16 to the data of L*'a*'b* f color 
space of a color processing input to check checking the precision of color processing It means 
investigating what has happened to the data of a color space. Therefore, change the data of an L*'a*'b* ? 
color space of arbitration to check on the three-dimension table 111, and it asks for the data of a 
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Y'M'CK' color space. The value is changed into the data of a YMCK color space using the 1- 
dimensional tables 1 12-1 15. L* a* b* outputted from a color printer 1 16 by making the data of a YMCK 
color space into a forecast-ed What is necessary is to predict the data of a color space and just to search 
for the color difference with the inputted data of L^a*!)* 1 color space. 

[0128] Next, the case where a forecast-ed is 1 -dimensional data is again explained as an example which 
applied the gestalt of operation of the 3rd of this invention using drawing 7 and drawing 8 . As 
mentioned above, drawin g 7 shows the color processing in the case of a color scanner input unit, and 
explained previously the example which decides the coefficients of the matrix masking section 105 to be 
the 1 -dimensional tables 102-104. Here, these coefficients are already decided once, color transfer 
characteristics change with the light source of a color scanner input unit, time fluctuation of CCD, etc., 
and it explains supposing the case where only the coefficient of the 1 -dimensional tables 102-104 is 
decided on again again in simple. 

[0129] Drawin g 8 shows the color processing coefficient decision method when not using a data 
precision function, and also explains a procedure here according to the procedure of drawing 8 
supposing the case where a data precision function is not used. The gray-scale patch which has 
measured lightness beforehand is made to read at the I/O live-data pair preparation production process 
of S81 first, the signal in the RGB color space of an output is measured, and it is lightness L* as live 
data of an input signal. A live-data pair with the data in the RGB color space as live data of an output 
signal is prepared. Although it is dependent on the precision of the prediction for which it asks, 10 to 
about several ten are usually suitable for the gap of gray scale, and when it also determines the 
coefficient of the matrix masking section 105, it can reduce the number sharply. Since the noise by 
causes, such as heterogeneity within the field of a color scanner and time instability, has ridden on the 
output signal at the time of reading, the device from which these noises are easy to be removed by the 
smoothing feature which the color property prediction method has, the device also in consideration of 
the heterogeneity of the gray-scale patch [ itself], etc. are the same as that of an above-mentioned 
example. 

[0130] Next, at the weighting function configuration / parameter setup production process of S82, the 
function configuration and parameter of weighting at the time of color transfer-characteristics prediction 
are set up. In this case, what is necessary is just to set up the weighting function equivalent to (8) types 
like (40) types, since an input is lightness. Here, the input signal of a forecast-ed and (Li) are the 
constants (Lj) decides the live data of an input signal, and the constant of standardization of (L0) of 
input signal space, and p decides function configurations to be. 
[Equation 11] 

W Uj = 1/((((I« - Lj)/L Q YY + 1) (40) 

[0131] It is as having mentioned above that it is what the parameter of weighting means (L0) and two 
kinds of constants of p, and adjusts the degree of a smoothing feature. In this case, since it assumes 
having determined the coefficient of the 1 -dimensional tables 102-104 and the matrix masking section 
105 once, the parameter of these weighting refers to the parameter then used, and should just determine 
it. 

[0132] Next, at the forecast-ed criteria production process of S83, two or more forecasts-ed for deciding 
the coefficient of the 1-dimensional tables 102-104 are set up. (In this case, gray scale, i.e., lightness L*, 
What is changing equally from zero to 100 is set up with a forecast-ed.) It will be satisfactory, if width 
of face of the step is made about into 100/256 when using a 8-bit table. 

[0133] Next, at the forecast calculation production process of S84, prediction is performed about two or 
more forecasts-ed set up at the front production process. The method will be entire identitas if the point 
that processings from S51 to S56 explained in drawin g 4 differ from the number of dimension of an 
input is removed. Therefore, explanation is omitted here. The forecast in the RGB color space to the 
gray scale which is a forecast-ed is obtained at this forecast calculation production process. 
[0134] Next, at the data election production process for color processing coefficient decision of S85, a 
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data pair required in order to decide the coefficient of color processing is elected from the data pair of 
two or more forecasts-ed and a forecast obtained even at the front production process. The criteria of 
election are whether the close forecast in a RGB color space is in the range treated with a color 
processor. That is, although the data of a RGB color space can be expressed with the unit of a reflection 
factor, the data with which the forecast in the RGB color space of a gray scale exceeds less than 0 and 
100% may usually be contained. What is necessary is to cut those data and just to elect the data pair to 
which the close forecast in a RGB color space is in 0 to 100% of range, since those data is unnecessary 
as data for color processing coefficient decision. 

[0135] Next, the coefficient of an actual 1 -dimensional table is determined at the color processing 
coefficient decision production process of S86. The 1st purpose of the 1-dimensional tables 102-104 is 
data of the RGB color space usually proportional to a reflection factor at the following matrix masking 
section 105 L* a* b* When changing into the data of a color space, it is changing into the data of the 
R f G'B' color space which makes an unit another unit which is not a reflection factor, for example, 
concentration, lightness, etc., so that conversion precision may become good. Moreover, when gray is 
inputted into a scanner, the 2nd purpose is changing so that the data of a R'G'B' color space may become 
the same value, in order to amend the difference in the sensitivity of the RGB sensor of a scanner. Here, 
the case where it changes into the unit of lightness is assumed. The coefficient of R, G, and the 1- 
dimensional tables 102-104 corresponding to B of each is the value of the forecast of R, or G or B to an 
axis of ordinate in a horizontal axis Lightness L* A value is taken, the data pair of a gray scale is 
plotted, the transform function from a RGB color space to a R'G'B 1 color space can be defined by 
approximating with the polygonal line between plots, and it can be determined by quantizing this 
transform function. 

[0136] At the last, the precision of color processing is checked at the color processing coefficient check 
production process of S87. Lightness L* of the scanner input of the gray scale which is color processing 
to check checking It receives and is L* a* b* after color processing. It means investigating what has 
happened to the output data of a color space. Therefore, gray scale L* of arbitration to check It considers 
as a forecast-ed, the output data of the RGB color space of a scanner are predicted, it is changed into the 
data of a R'G'B' color space using the 1-dimensional tables 102-104, the matrix masking section 105 is 
used further, and it is L* a* b*. It changes into the data of a color space and is lightness L*. What is 
necessary is just to search for a difference. 

[0137] As mentioned above, if the color transfer-characteristics prediction method of this invention is 
used and 1-dimensional data processing is performed also when deciding on only the coefficient of a 1- 
dimensional table again again, realizing using few live-data pairs is possible. 
[0138] Although this example explained the case where a data precision function was not used, of 
course, a processing process as shown in drawing 10 can also perform redeciding the coefficient of a 1- 
dimensional table using a data precision function. 

[0139] In addition, coefficient Wlij of weighting which made the input signal the forecast-ed and was 
shown by (40) formulas by making an output signal into a forecast here Although used prediction 
processing was performed An output signal is made into a forecast-ed and an input signal is made into a 
forecast-ed. A part of output signal and input signal Conversely, a forecast-ed, If the above-mentioned 
successive approximation method is used when performing prediction processing by making a part of 
remaining input signals into a forecast (for example, a case like the color printer shown in drawing 9 ), it 
will be Wlij as a coefficient of weighting. It can use and prediction processing can also be performed. 
Of course, in addition to it, a data precision function may be used. 

[0140] Thus, in various cases, the color processing coefficient decision method of this invention can be 
applied, when determining the coefficient of the color processing about a color picture input unit or a 
color picture output unit, or when checking the precision of the coefficient of the color processing for 
which it opted. 
[0141] 

[Effect of the Invention] According to this invention, the live data of I/O of a color picture I/O device 
can be processed statistically, and the input signal over prediction of the output signal over prediction of 
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color transfer characteristics, i.e., the input signal of arbitration, and the output signal of arbitration can 
be predicted, or a part of remaining input signals can be predicted from a part of the output signal and 
input signal of arbitration so that clearly from the above explanation. By this, color transfer 
characteristics can be predicted without being dependent on the class of color picture I/O device, and the 
time and effort which investigates and models the physical property of equipment can be saved, and it 
can predict easily only by deciding the parameter of weighting to any color picture I/O devices. 
Moreover, even if the noise is contained in the live-data pair of the I/O to be used, by determining the 
parameter of weighting with the degree of the noise, high predictability can be obtained and the 
coefficient of color processing of a color picture I/O device can be decided correctly. Furthermore, since 
the prediction method is easy, it can predict at a high speed and there are various effects ~ the 
coefficient of color processing of a color picture I/O device can be determined quickly. 



[Translation done.] 
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